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Summary Temperature and humidity measurements of air in and above a drying fluidized mass 
of silica gel are given. It is proved that perpendicular to the flow of air the temperature and the 
humidity show no gradient ; that the temperature of the silica gel is practicaliy equal to the 
temperature of the leaving air, being independent of the place in the apparatus ; that mass and 
heat transfer takes place in a very small layer only and finally that heat is transferred to the 
solid substance by the supporting sieves as well as by the air directly. 


Résumé—Les auteurs ont effectué des mesures de la température et de humidité de l’air dans 
et au-dessus d'une masse de silicagel séchéante et fluidifiée. Ils prouvent que, dans un sens 
perpendiculaire au courant d’air, la température et humidité ne présentent pas de gradient ; 
que la température du silicagel est pratiquement égale a la température de lair sortant puisque 
cette température est indépendante de l'endroit dans l'appareil ; que la transmission de la 
substance et de la chaleur ne s'’effectue que dans une couche trés mince et enfin que la chaleur 
est transmise & la substance solide tant par les cribles supportants que directement par lair. 


INTRODUCTION 
In the course of our studies on the drying of a 
fluidized mass with the aid of air it proved to be 
of importance to measure the temperature and 
the humidity of the air in the fluidized mass, 
as well as the surface temperature of the solid 
substance. As far as we could see neither the 
measurements of KeTrenrING et al. [1] nor those 
of Leva et al. [2] gave the correct air conditions 
in the bed. We therefore tried a different 
technique in obtaining the desired data. The 
results of our experiments will be given below. 
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MetHop oF INVESTIGATION 

A closely graded quantity of silica gel is fluidized 
by means of a current of air of known quantity, 
temperature and humidity in a_ cylindrical 
container, — with one or more sieveplates - 
insulated against heat flow by means of a jacket. 
To measure the conditions of the air in the silica 
gel a tube of rather small diameter, adjustable 
in a horizontal and vertical direction, was inserted 
into the fluidized mass. Through this tube air] pig 1. Sketch of the container of the fluidization 
could be sucked by means of a waterjet pump apparatus. 
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at an adjustable velocity. The tube could be 
sealed by means of sicves of different materials 
to prevent the silica gel entering the small tube. 
The following measurements were made: the 
temperature and the humidity of the incoming 
air and of the air leaving the fluidized mass ; 
the temperature and humidity of the air in the 
inner tube (on a short distance above the sealing 
sieve) and the temperature in the fluidized mass, 
measured in the same plane 
to the direction of flow 
the inner tube. 


perpendicular 
as the temperature in 





t, = temp. of air leaving the fluidized mass 
measured some cm above the surface of 


this mass. 
t,, = temp. in the fluidized mass. 
t, = temp. of the air in the inner tube measured 


just above the small sealing sieve. 

The humidity was measured with a Foxboro 
dewcell recording meter and expressed as the 
dewpoint. 

The humidity of the incoming air was 
always kept constant (dewpoint in different runs 
ranging from 10 to 13°F (+ 0-2°F)) by passing 


























A sketch of the fluidization apparatus used, 
drawn as far as possible to scale, is given in Fig. 1. 
The cylindrical vessel and jacket were made of 
Pleviglass, the inner tube of glass. The main 
sieve used was of celluloid. Temperatures were 
measured with thermocouples of copper-constan- 
tan with a small welded point and as great a 
length of wire as possible in an equithermal 
plane, via the compensation method. An excep- 
tion was made for the incoming air, the tempera- 
ture of which was measured with a smal! mercury 
thermometer. The temperatures measured were : 
t, = temp. incoming air, measured some cm 








below the sieves. 


dry air (tapped trom the compressed air system, 
after the removal of liquid water) through a bed 
of blue silica gel. The temperature of this air 
could be adjusted by means of an electrically 
heated wire with adjustable current. The 
quantity of air was measured with a flowrator. 

The humidity of the air sucked through the 
inner tube and of the leaving air were again 
measured with the Foxboro meter the quantity 
of the air passing the inner tube measured by 
means of a capillary tube flow meter, later with 
a flowrator. 

A sketch of the whole set-up is given in 
Fig. 2. A is the flowrator for the incoming 
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air, B and G open manometers, C and M are 
regulators to adjust the quantity of air, D 
contains the blue silica gel, E is the heater with 
wire, resistances and ammeter F, H is the 
jacketed fluidization apparatus of a somewhat 
different kind as given in Fig. 1. (no changeable 
sieves). J is the inner tube ; at J a tube could be 
attached to measure the humidity. K contains 
the dewcell of the humidity meter; L is a flow- 
meter ; O are two india-rubber tubes and finally 
at N the vacuum pump was connected. E and 
H were insulated by means of cottonpads. 

As variables in the measurements were 
chosen, time, the places of measurements in 
and above the fluidized silica gel, the flow rate 
and temperature of the incoming air, type and 
number of sieves used to support the silica gel 
and of the small sealing sieve(s) for the inner 
tube. 

Silica gel of about one size, namely 0-3- 
0-5 mm, was used with a constant initial water 
content of 42-8°%, of water calculated on the dry 
substance. Reproducible initial amounts of 
water in the gel could be easily obtained by 
fluidizing the gel in a separate apparatus with 
moist air, humidified in two packed towers 
filled with Raschig rings and water. After the 
towers, the air passed a separator to catch the 
entrainment. The whole set was placed in a 
thermostat at 20°C. A sharp indication for 
the endpoint of the humification was a constant 
temperature in the gel. 


Tue Resutts Osratnep 

(a). The entrance effect (e) 

A typical set of curves of a normal run is given 
in Fig. 3 where t, and ¢, are given as a function 
of the drying time. The temperature of the 
leaving air initially falls off very rapidly and then 
increases again rather rapidly, becomes constant 
during a certain period (which will be called the 
constant rate period of drying throughout this 
paper) and finally rises again reaching assymtoti- 
cally the temperature of the incoming air. In 
Fig. 3 is also included the temperature in the 
inner tube placed in the centre of the vessel 
(as in most measurements) without the sealing 
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sieve, directly on the apparatus sieves; ¢; at 


0mm above the sieves (= t,,.). The amount of 
air sucked through the inner tube was equal to 
(4, 
(Di)? 
dia meter of the inner tube, D, the inside diameter 
of the apparatus. 


total amount of air, d, being the outside 
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Fig. 3. t,, t; and | as a function of the time of drying. 


It is clear from Fig. 3 that t,, <<t,. From 
this it can be concluded that the heat transfered 
from the incoming air to the fluidized substance 
is transported in at least two ways viz. (a) 
directly and (5) via the supporting sieves by contact 
between substance and sieves (see also the 
measurements of Leva et al. [2]). In our case 
even with an adiabatic apparatus there will 
be a heat exchange as well with the walls of the 
apparatus, changing in sign during the operation. 


Neglecting these and other heat exchanges 


the ratioe = ‘ — is a measure for the transfer 
of heat to the gal by means of the supporting 
sieves. We found that, within the accuracy of 
our determinations (0-1°C) and excluding the 
first 5 minutes of every experiment for obvious 
reasons, € is independent of the time of drying. 
As an example in Table I the results for three 
experiments differing in quantity of air are 
given. 

It is clear that e besides being a function 
of the air velocity as shown in Table 1, also 
must depend on the material of which the sieves 
are made, the thickness of the sieves and the 
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Table 1. ¢ as a function of the time of drying. 
: | 
time (min) > 5 10 20 30 40 50 | 60 «)=siés7 sO 
" Viux 36-61/min | 0-48 | 0-47 | 0-45 | 0-48 | 0-48 | 0-48 | 0-48 | 0-48 | 0-468 
' . ' ee | — 

: Flux 64-2 1/min 0-28 | 0-29 | 026 | 0-25 | 0-25 | 0-26 
Flux 84-1 1/min | o-20 | 0-21 | 022 | 0-20 0-19 018 om 
time (min) > 0 100 | 110 | 120 | 130 | 140 | 150 | 160 | 170 
Flux 36-4 1/min 043 O42 | 0-44 0-43 0-48 0-41 0-42 | 0-42 | 0-40 

e — -_ = 


Flux 64-2 1/min 0-25 029 | 0-27 


Flux 84-1 1/min 0.27 


0-24 | 0-28 | 0-27 


‘. meme 6 —_}—__ 


| 





The temperature of the incoming air in all cases was about 20°C, the initial amount 
of wet silica gel 95 g. Drying was in all cases continued until a water content of about 10%. 


free passage of the sieves. This holds good for 
the supporting sieves as well as for eventual 
sealing sieves of the inner tube. In Fig. 4 some 
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Fig. 4. The entrance effect (¢) as a function of the 
amount of air and with different systems of sieves. 


data obtained with experiments using differing 
systems of sieves have been collected. e was 
measured at the “ constant rate " period, because 
this facilitates the measurements. In all cases 
the inner tube rested on the supporting sieves in 





the centre of the apparatus. In the direction 
of flow the air passed respectively : 


apparatus inner tube 


(a) Nickel gauze /celluloid sieve / — 


_— 


(b) copper sieve /celluloid sieve /nickel gauze — 


(c) nickel gauze 
‘celluloid sieve / nickel gauze 


(d) celluloid sieve /copper sieve celluloid sieve 


(e) copper sieve /celluloid sieve 
/copper sieve celluloid sieve 
(f) copper sieve /celluloid sieve 
/nickel gauze celluloid sieve 
(g) copper sieve /celluloid sieve 


/nickel gauze celluloid sieve 


The celluloid sieve of the apparatus was 
always the same, 2mm thick with 33 holes 
of 0-4mm diameter per cm?*, the small celluloid 
sieves applied to the inner tube differed in free 
passage — which was rather small — decreasing 
from d +g. The nickel gauzes of 0-065 mm 
thickness had a rather large free passage, con- 
taining 1800 holes of 70, diameter per cm, 
whereas the copper sieve was fabricated by 
drilling 80 holes of 0.5mm diameter per cm? 
in a copper plate of 0-5mm thickness. No 
further quantitative data about the sieves used 
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will be given here because no effort to calculate 
the entrance effect has been done. It can be 
mentioned however that under simplifying 
assumptions regarding the type of flow under 
and through the sieves, it has been possible to 
calculate the decrease in temperature over the 
sieves, assuming heat is transfered in the sieves 
by conduction. This calculated decrease was of 
the same order of magnitude as found. 





fluidized systems, which otherwise looks rather 
inviting in view of the homogenity as far as the 
dimensions of the apparatus is concerned. 


(b) Temperature and humidity measurements in 
planes perpendicular to the direction of flow. 


It seemed interesting to know whether there 

would be a temperature or a humidity gradient 

in planes perpendicular to the direction of flow 
i.e. in horizontal planes. 
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‘as Above the fluidized mass the 
temperature as well as the humidity 
of the air shows no gradient in a 
horizontal plane as far as the possi- 
bilities of measurement could show, 
as can be seen from the Figs. 5 and 
6. In Fig. 5 two experiments are 
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given, I with 95g of wet gel and 
the temperature measured 2 cm 
above the fluidized mass at three 
places as indicated; II with 50g 
of wet gel, measured at three places 
4cm above the mass. In both cases 
the temperatures plotted against 
time, show that there exist no differ- 
ences between the three places. 
In Fig. 6 the dewpoints 
measured in the centre of the 
apparatus (0) and 6mm from the 
wall (z) as well as the temperature 
t, are plotted against time and again 
show that there is no gradient. In 








The temperature of the air above the gel measured 


at different places in a horizontal plane. 


The data collected in Fig. 4 prove that 
indeed e depends on air velocity as well as on 
the material, thickness and free passage of the 
sieves. The entrance effect will be small with 
very thin sieves with a large free passage and of 
a material of low conductivity. Also the higher 
the air velocity, the smaller will be the entrance 
effect. Within the range of velocities used e 
x (velocity)°* 

The entrance effect strongly hinders the 
study of heat and mass transfer problems in 





this last case 95 g of wet gel was 
used, the plane was 2cm above 
the fluidized mass. In all the three 
experiments about 65 litres of air 
passed through per minute. 

In accordance with the above results the 
air in the gel also showed no appreciable gradient 
in a horizontal plane. 

In Table 2 some of the results obtained are 
given. As a measure for the air temperature is 
given e,, expressed in % based on ¢; a distance 
zmm from the supporting sieves. At the same 
time ¢,, (the temperature in the fluidized mass) 
was measured and arbitrarily will be represented 
t—t 


, > _t—t,, 
in the number ¢,,, in %, based one, = *— 


te —t 
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Fig. 6. The dewpoint of the air above the gel measured 
at different places in a horizontal plane. 


100%, in which ¢,, is the temperature in the 
mass at a distance of z mm from the supporting 
sieve. All numbers were based on measurements 
during the constant rate period. e, and e¢,., 
are given as a function of the distance from the 
centre of the apparatus expressed in parts of the 
radius of the apparatus. 


Table 2. ¢;, and ¢,,, as a function of the distance 
from the centre of the container. 





Distance 
in radii 


Distance 
in radii 


Distance 
in radii 


. ; o 
finy im % 


Distance 


in radii | 035 O35 O54 0-93 


my in % | 82 77 : | 95 





of the 


container. 





All experiments were carried out at the same 
velocity of 64 1/min, 95 g of wet gel was used, 
the temperature ¢, = 20, 5°C, the apparatus 
had a nickel gauze under the celluloid sieve, 
the inner tube was sealed with a nickel gauze. 

In two other experiments using 95g and 
50g of wet silica gel, using the nickel gauze 
celluloid sieve and the inner tube resting on the 
celluloid sieve, at the distances of 0 and 0-54 
radii from the centre was found ;: e,, resp. 24% 
and 24%; ¢,,,: 81% and 78%. From these 
examples it is clear that the air temperature in 
the fluidized mass is also constant in a horizontal 
plane. The temperature ¢,, as is evident, is 
lower than the air temperature in the gel, moreover 
It can be concluded therefore 
that this temperature (as has been done by others) 
cannot be taken for the temperature of the air 
in the gel, because it represents a mean tempera- 
ture of air and surface temperature of the gel ; 
moreover the surface temperature of the silica 
gel must be lower than the air temperature and 
finally this mean temperature will depend on the 
frequency of contact of the thermo-couple with 
the gel. There is a tendency that the mean 
temperature at the wall is lower than in the 
midst of the apparatus. This would mean a 
higher concentration of gel at the wall, which 
again would tally with the observed movements 
in the fluidized bed, the particles going up in the 
centre and coming down near the wall. The 
values of t, however are not very reliable, 
and as the authors see it, ought not to be used. 

It should be remarked here that the authors 
have found a definite difference in pressure 
between the air in the midst of the apparatus and 
at the wall, measured in a horizontal plane 
1 cm above the sieve, the air in the midst showing 
a smaller pressure. Differences in pressure- 
difference with the air outside the apparatus 
up to 8% could be found. 

It can be concluded from this paragraph 
that within the accuracy of the determinations 
the air in the apparatus in and outside the 
fluidized mass shows no temperature, nor a 
humidity gradient in a horizontal plane. The air 
quite near the wall was not included in the 
observations, therefore also not in the conclusion. 


it is not constant. 
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(c) The temperature of the fluidized gel. 

As has been said in paragraph (b) the mean 
temperature in the fluidized mass is lower than 
the air temperature ¢; in the gel and therefore 
the surface temperature of the silica gel must 
be lower than this temperature t,. 

In the next paragraph (d) it will be discussed 
why the authors are of the opinion that the 
temperature ¢,, measured in all cases while a 
quantity of air is sucked through the inner 


(d,)? 


. x the total quantity of air 
(Dj)? 


gives a fair approximation to the air temperature, 


tube equal to 





thermocouple was moved into the bed to a 
height of about 2cm above the sieves, the air 
supply stopped and after about 10 seconds the 
temperature read. Prolonging the time to about 
one minute did not change this temperature in a 
measurable way. 

Then the couple was moved upward, the 
air supply opened again and the measurements 
repeated. Two typical temperature-time curves 
are given in Fig. 7. The temperature t, (above 
the bed) is given as points (-), the temperature 
in the bed by points (0). These points indicate 
that the temperature of the gel is almost equal 
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Fig. 7. 


—--— = 


The temperature of the gel (0) compared with 


the temperature ¢, of the air above the gel (.). 


if a thin nickel-gauze of large free passage is used as 
a sealing sieve. (d, = outside diameter of the inner 
tube, D; = inside diameter of the container). 

In this paragraph it will be discussed how 
the authors have tried to evaluate the surface 
temperature of the fluidized mass. 

During a normal drying experiment the 
inner tube was removed. A thermocouple which 
could be replaced in a vertical direction was 
placed in the axis of the apparatus. The tem- 
perature ¢, was measured above the fluidized 
mass as a function of the time. At intervals the 





to the air temperature above the bed, with a 
possible difference of the accuracy of the tempera- 
ture measurements, viz. 0-1°C. The small 
amount of stagnant air present has practically 
no influence on the temperature. 

Near the supporting sieves the procedure 
described above could not be followed, because 
the hotter sieves transported heat to the ther- 
mocouple too quickly, hindering the measure- 
ments. By following the behaviour of a ther- 
mocouple inserted into the fluidized mass indicat- 
ing t,, at a height of 14mm above the sieves, 
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it was however observed that after closing the 
air supply the temperature of ¢,, fell rapidly 
and then rose rapidly again. The lowest tempera- 
ture observed was only some tenth of a degree 
centigrade higher than t, above the gel. 

Because of the very small particles used, 
causing an equalising time of the temperatures 
in the gel to differences < than 0-1°C of some 
tenth of a second if initially these temperature 
differences are of the order of 10°C, as can be 
calculated with the aid of the Fourier number, 
it will be clear that the temperature of the gel 
measured is a mean mass temperature and not 
a surface temperature. 

On the other hand for the same reason the 
temperature differences in the gel will be small, 
and with these small particles therefore a distince- 
tion between surface temperature and mean mass 
temperature real practical meaning. 
Therefore the temperature measured is called the 
temperature of the gel. 

From the observations the conclusion seems 
warranted that under the used in 
our experiments, after a short time of drying 
the temperature of the fluidized gel at any time 
is almost equal to the temperature of the air 
leaving the bed and independent of the place 
in the bed. 

The actual temperature 
between the sieves and the gel in contact with the 
sieves could not be measured. 


has no 


conditions 


drop _ existing 


(d) The temperature and the humidity of the air 
in the fluidized bed in a vertical direction 

By moving the inner tube in a vertical direction 
the temperature and the humidity of the air 
in the fluidized mass at different heights could be 
measured. The measurements were all carried 
out with the tube in the axis of the apparatus 
and to facilitate the measurements during the 
constant rate period of drying. 

As said before it still has to be proved that 
the air in the tube in reality gives the air con- 
ditions in the fluidized mass as wanted. To 
prevent the removal of solid substances the 
inner tube had to be sealed with a sieve. For 
this we chose a very thin nickel gauze of large 
free passage (see section (a), p. 99). 





To prove that the temperature and humidity 
measured closely approximate the desired values 
the following sets of experiments were carried out. 
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Fig. 8. ej, of the air sucked through the inner tube, 


(placed 1mm above the copper sieve) as a function of 
the amount of air. 
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Fig. 9. ¢jg and the dewpoint of the air sucked through 
the inner tube (placed 3mm above the celluloid sieve) 
as a function of the amount of air. 


In the first set the humidity and t; were measured 
with different amounts of air sucked through 
the inner tube. The results obtained are collected 


in Figs. 8 and 9 in which ¢ = so and the 
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dewpoint are given as a function of the amount 
of air passing the inner tube, expressed in litres 
per minute and in °%, of the total amount. 
Measurements were again carried out at the 
constant rate period. The tube was 1mm 
(Fig. 8) and 3mm (Fig. 9) respectively above 
the sieves. 

Both the curves obtained show a maximum. 
The amount of air passing the inner tube normally 


; ; d.\? 
used in our experiments was as stated ae x 
| 
100%. As d, 1-26cm and D; = 6-0 cm, this 
amounts to 44%. Both maxima lie approxi- 
mately at this velocity (see dotted vertical lines) 
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Fig. 10. 4; as a function of the reciprocal number of 
holes (1/n) in the sealing sieve, meusured at different 
heights. 


At velocities higher than 44% part of the 
substance is sucked against the nickel gauze 
causing less movement of this part, resulting 
in an accelerated drying giving relatively higher 
temperatures and lower humidities. At lower 


Table 3. 





has been shown this sieve can have influence 
on the temperature, by the substance 
moving along its surface and taking up heat, 
decreasing the temperature of the sieve and 
therefore of the air that has passed through it. 

A set of experiments was now carried out, 


solid 


using sieves of celluloid of the same thickness 
(0-25mm) and with an increasing number of 
holes of the same diameter (0-2mm). The 
temperature ¢, obtained under the same con- 
ditions, at an amount of 4.4% from the total 
of air passing the inner tube, and measured at 
the constant rate period, is plotted in Fig. 10 
against the reciprocal number 1/n of the holes. 
The values of t; extrapolated from this figure 
to the maximum of holes possible, calculated 
from the total area of the tube (see dotted 
vertical line), are given in Table 3 for four 
different distances from the supporting sieve. 
In the same table t; as measured with the normal 
nickel gauze as a sieve and for the same conditions 
are given. 

From these data it can be concluded, that 
the nickel gauze used does not decrease the 
temperature of the air if the tube is placed in 
the gel. 

It must be remarked however that if the 
distance of the inner tube to the supporting 
sieve = 0mm, so that the nickel gauze is in 
contact with this sieve, a decrease in the tempera- 
ture ¢; is observed, which amounts to about 
0-3°C under the conditions described. 

From these two sets of experiments it is 
clear that ¢; as measured in fair approximation 
gives the temperature of the air as desired. 

It may be remarked here, that inevitable 


Determination of t; for different sieves. 





Distance from supporting sieve in mm 4 2 + 6 


/; extr. for total free tube in °C 


1, determined with nickel sieve in °C 


15°8 
15°8 


16-9 
17-0 


216 20-1 
215 203 





velocities than 4-4°%, the air will have had a 
less intensive contact with the fluidized mass 
resulting again in relatively higher temperatures 
and lower humidities. 

In the second set of experiments the influence 
of the nickel sieve on ¢; was measured. As 





disturbances in the bed caused by the insertion 
of the inner tube, seem to be of no great influence, 
if the distance of the tube to the supporting 
sieves is not too small. 

The results of measurements of the tempera- 
ture and the humidity of the air in a vertical 
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direction in the fluidized mass for the case in 
the amount the 
amount of wet silica gel was 95 g and the celluloid 


which of air was 65 1/min, 
sieve was used, are represented in Fig. 11, in 


which ¢,, t,, and ¢, and the dewpoint are plotted 





bed, whereas the rest of the bed is inactive 
in this respect. For drying apparatus this is 
important as far as the height of the drying 
the of the 
concerned, 


bed and _ therefore construction 


apparatus is 
That 


the dewpoint 
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after 6mm has not yet 
fully reached the maxi- 
as es mum value may be due 
' either to some lag in the 


; 


measuring device, and, or 
to the diffusion of water 
vapour from the pores of 
the solid substance to the 
main stream of air. 

The differences in the 
form of the t; — distance 


| , above ge/ | 
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’ pitial (10 °F) 
miansh 
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lines in both cases may be 


Dewpoint ——e 


caused by the difference 
in sieves used, causing a 
different movement of the 
particles over the sieve 





t surface. 








Fig. 11. 


The temperature and the humidity of the air as a 


1g 18 . ° 
The conclusion given 

is supported by experi- 

ments done with different 


mm 


function 


of the distance above the supporting celluloid sieve. 


the distance from the sieve in mm. 


The inner tube was in all cases sealed with a 


against 


nickel gauze. The total height of fluidized mass 
about 60mm. Point A in Fig. 11 gives 
the temperature ¢; if the inner tube is not sealed 
with a nickel gauze (see above). In Fig. 11 also 
the temperature ¢, and the dewpoint of the 
incoming air, as well as the dewpoint of the 
leaving air are given. 


was 


In Fig. 12 the temperature measurements 
obtained with the copper sieve mentioned before, 
instead of the celluloid sieve are represented. 
The conditions are the same as regards quantity 
of air and amount of substance. The temperature 
given in Fig. 12 is a temperature measured with 
a thermocouple pressed under and against the 
sieve. A same result was obtained in the experi- 
ments collected in Fig. 11 (not given). 

The of these experiments is 
that heat exchange and mass transport take 
place practically only in the first 4-6 mm of the 


conclusion 





amounts of substance, 





=) 
° 























am 7 = 

REZ iA eg & AK, 

Distance above sieve mm 
—_—_ o& 





Fig. 12. The temperature of the air as a function of the 
distance above the supporting copper sieve. 
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dried under the same conditions (amount and 
temperature of air, and sieves used). The results 
are collected in Fig. 13 giving the temperature 
t, (x), the temperature ¢,, (+-) and ¢; (-) measured 
0mm above the sieve, and the dewpoint (0) 
of the leaving air as a function of the amount 
of gel, all measured during the constant rate 
period. Because the measurements at low quanti- 
ties of substance are extremely difficult and not 
very reliable, due to the short period of constant 
drying conditions, relatively bigger influences of 
the walls (causing an increase in t, while ¢, is 
constant) and poor fluidization (see ¢,,), Fig. 13 
warrants that heat transfer 
takes place for the greater part only at the 
surface of and in a rather thin layer above the 
supporting sieves. This layer was in our experi- 
ments about 5mm. thick. 


the conclusion 
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Fig. 13. t, ty, t,, and the dewpoint of the air as a function 


of the amount of gel. 


As far as heat transfer alone is concerned, 
it can be calculated from data given in the 
literature that under the conditions mentioned, 
the heat exchanging layer is about 2-3 mm 
thick if a difference of 0-1°C is assumed between 
leaving air temperature and surface temperature. 


CONCLUSIONS 
Within the limitations given by the conditions 
used the following conclusions can be given for 
the drying of a fluidized bed of silica gel with 
air in an adiabatic apparatus using a supporting 
sieve and vertical flow of air : 





(a) The air temperature and humidity are 
independent of the place in the apparatus 
in a horizontal direction ; 

(b) The temperature of the gel is independent 
of the place in the apparatus and equal to 
the temperature of the air leaving the bed ; 
(ce) Heat is transferred to the gel in at least 
two ways; by means of contact of the 
fluidized substance with the sieves and by 
means of contact with the air. The ratio 
of the amount of heat transferred by the 
sieves over the total amount transported 
is independent of the time of drying. The 
ratio decreases with increasing air velocity 
but increases with decreasing free passage 
in the sieves and with increasing thickness 
of the sieves. 

Practically all heat and mass transfer 

takes place in a rather thin layer above 

the sieves. The depth of this layer does not 
depend appreciably the amount of 
substance fluidized. 


on 


NOTATION 


la 


wa lio ' : ' 
———  . (No dimension). 
— ty 


¢ = The entrance effect = 


la—liz } 
— (No dimension) A 


- certain ratio of 
tg—ty 


temperature difference 
t measured z mm above 
; the sieves. e.g. €;), €jg- 


= (No dimension). 
fg — % 


Outside diameter of inner tube (cm). 
= Inside diameter of the container (cm). 
= Temperature of the air entering the apparatus 
(°C). 
t, = Temperature of the air leaving the fluidized mass 
(°C). 
= Mean temperature of the fluidized mass and air 
measured z mm above the sieves. (°C). 
tj, = Temperature of the air in the inner tube of the 
apparatus, measured z mm above the apparatus 
sieves. (°C). 
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Hydroextraction V: Radial distribution of pressure in homogeneous 
and banded cakes 


M. M. Harun: and J. ANDERSON StTorKRow 


Chemical Engineering Laboratory, Applied Chemistry Department, College of Technology, Manchester 
( Received 7 July 1953) 


Summary —Hydroextraction and filtration permeabilities have been measured in a 9 in. dia. perforated 
basket centrifuge for cakes formed from layers of kieselguhr of high, and of starch and of chalk of low 
permeability. Conditions that lead to discontinuity in the development of centrifugal head in the liquid 
within a heterogeneous cake are absent from cakes of one material. The results justify the assumption 
that the latter cakes run with their void space full of liquid and that atmospheric pressure is transmitted 
to the outer surface of the cake in a hydroextractor. This flow pattern is the fundamental basis of the 
proposed equations for hydroextraction. 


Résumé— Dans un panier d'essoreux centrifuge de 225 mm. de diam. T'auteur détermine les perméabilités 
pour l'essorage ct pour la filtration de couches alternées de kieselguhr (haute perméabilité) et d'amidon ou 
de craie (faible perméabilité). Avee des gateaux homogtnes, on ne retrouve pas les conditions qui con- 
duisent & létablissement d'une discontinuité dans la force centrifuge agissant sur le liquide présent dans 


un gateau hétérogéne. 


Ces resultats justifient Phypothése qu'un tourteau homogéne tourne avec ses pores remplis de liquide 
et que la pression atmosphérique progresse vers la surface externe dans l'essoreux. Cetter torme d écoule- 
ment est prise comme base pour les equations de lessorage, proposées por l'autcur. 


It has been shown elsewhere [1], [3] that the flow of 

a liquid through a hydroextractor cake of constant 

permeability K follows the relationship 
4n3°n? XK (r,? —r,*) 

" we " loge To /Te 
when the centrifugal head generated in the system is 
resisted only by viscous shearing forces in the laminar 
motion of the fluid within the cake. 

It has been shown [2] that the permeability of the 
innermost layer may be 25% greater than that of the 
outermost layer of a thick cake. A cloth-cake inter- 
face resistance, R, was assumed to affect both processes 
of hydroextraction and filtration. Thus the general 
form of equation (1) for a cake composed of layers 


becomes : 


(1) 


, 
log, -° 
> 


4n° nn? X (7? — r,?) : 


#E&-q K, 


+ 


These equations were based on the equality between 





the centrifugal head, P., and the energy dissipated in 
friction through the cake, P,, assuming that the 
pressure is atmospheric at r, and r, and that the 
cloth support is running “ empty,” i.e. without the 
void space in the cloth being completely filled with 
liquid [2]. 

It is of value to study the pressure distribution 
through a cake and the effect of a highly permeable 
layer in a cake, in order to assess the validity of the 
assumed flow mechanism. 

The pressure distribution will first be considered 
for a case conforming to equation (1), the modifica- 
tions appropriate to the system for equation (2) being 
considered later. 

The pressure terms are defined for a uniform cake 


(2a n)* 
sig is 


1 wq. r 
. ——z= log, 
K 227X "6 '. 


(r,* ae ry”) 


(3) 
P, = (4) 


At a constant speed of revolution, P., will be 
proportional to (r,? —r,*), and is zero at r =r;,. 
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Curves 4, B and C in Fig. 1 represent the values of 
Pat all points between r, and r, ; for r, = 0, 8, 8-8 
em respectively. P, is zero at r = r,, is proportional 
to g, and for a constant q follows logarithmic curves 
such as a, b and c. The curves a, b and c in Fig. 1 are 
for the value of g appropriate to each r, when P, = P, 
atr =r,. (P, — P,) represents the residual pressure, 
or manometric head, at any point on the fluid within 
the cake at r, << r < ro, and is zero at the cloth-cake 
interface, where the liquid is assumed at atmospheric 
pressure with the cloth running “ empty.” 























° 





Fig. 1. Variation of P, and Py, within a cake of uniform per- 
meability at n’ = 1,200 r.p.m., with different r,. 


In Fig. 2 (P, — P,) is plotted against r. From 
both figures it is seen that (P.— P,) is positive 
within the cake for r, = 0 or 8-0 cm, and is negative 
throughout the cake when r,; = r,, the liquid in the 
cake pores then being under a pressure lower than 
atmospheric. Fig. 3 represents this case on a large 
scale for a thick cake in the 9 in. diameter hydro- 
extractor at the upper limit of speed in the present 
investigations, which corresponds to the maximum 
centrifugal fields used in practice. The maximum 
suction head developed is 150cm water. As the 
distribution of actual K variations and the presence 
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RADIUS TF cm 
Variation of (P.-Pp) within a cake of uniform per- 
meability at n’ = 1,200 r.p.m., with different r,. 


Fig. 2. 


o. R tend to reduce this head, it seems likely that the 
maximum suction head developed in a hydroextractor 
cake in aqueous liquor is about 150 cm water. This 
manometric head is not low enough to cause boiling 





Nn 
N 





@ 
cm. water 





10 
RADIUS, F cm. 


, Pxio® 


PRESSURE 











10 
RADIUS, PF cm. 


Fig. 3. Variation of P, and Py within a cake of uniform per- 
meability at n’ = 3,000 r.p.m., 7, = Tr, = 8-46 em. 
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of water within the cake at atmospheric temperatures 
and hence there is no possibility for a vapour dis- 
continuity in the liquid which would invalidate the 
proposed flow equations. 

> , 
(P. Py) 


will be positive throughout the cake, but as r, 


When r, is considerably less than r,, 


increases, a lower than atmospheric is 


obtained at the outer part of the cake, as in Fig. 4, 


pressure 


extending progressively inwards as r,, tends to r,. 

If the P,, curve in Fig. 8 is recalculated taking 
into consideration the radial change in permeability 
[2] from inner to outer layers of the cake as 1-25: 1, 
the curve approximately coincides with P.. With a 
further resistance to flow, R, at the cloth-cake inter- 


the P, 


curve lies below that for P. at all r 
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Pp) within a cake of uniform per- 
3,000 r.p.m., with different r,. 





Fig. 4. Variation of (P, 


meability at n 


values between r, and r,. It is thus highly improbable 
in practice that any part of a hydroextractor cake 
It is thus 


improbable that either vapour or air bubbles will be 


will be under sub-atmospheric pressure. 


produced from the liquid within the cake to cause 
discontinuity in the liquid and hence in the pressure 
distributions. 

In the derivation of equation (1) it was assumed 
that the cloth layer behind the cake is not filled with 
liquid and hence does not affect the process of hydro- 
extraction. The earlier study of layer permeability 
[2] provided a technique for assessing the effect of 
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such a layer of highly permeable material forming 
part of a hydroextractor cake. 


EXPERIMENTAL TRCHNIQUE 


In a 9 in. dia, perforated basket centrifuge, one twill and one 
duck weave cloth layers were used as the support for the cakes. 
The cloths were spun wet at the highest speed before forming 
the cake. The cake material was fed to the basket as a thin 
aqueous slurry, the basket spinning at low speed. The surface 
of the cake was gently disturbed under water to help the for- 
mation of a vertical surface. This was followed by running the 
submerged cake for one minute at a high speed (2,000 r.p.m.). 
[4] Drainage time tests were done by measuring the time 
required for the liquid layer within the cake to move from 
rz, tO Tre at a specified basket speed (mainly 1,050 r.p.m. in 
this work). The probe arrangement [2] (Fig. 1) was used in 

(continued on the opposite page) 





Table 1. 
In hydroextractor basket : 
200 K + (100 K) to 700 K. 
3200S + (1008S) to 1.1008. 
200 K + (1008S) to 7008. 
3200S + (100 S) to 7008, 
3508 + 200 K 1008 
3505 + 100K + 2008 + 1508. 
3505 + 300 K + 2008 + 1508. 
2005 + 200 K + (1003S) to 5008. 


Layer Cakes. 


(w) 


+ (100 K) to 300 K, 
1008 150 S. 


In small spinning cell : 
9 2K + (2 K)to WK, 
0. 5C 4+ AC)two dC, 
11 25C4+4K 4+ 2-5C., 
12. 830 +4K + 2C, 
3. 40 +4K Bc. 
14. 5C +05K. 
5. 5C +2K. 
16. 5C + 4K, 
17. 2K+5C, 
18, 4K 4 5C, 
19% 6K 4 5C, 
20. 25C +G 
21. 4C4+4G4 
22. 5C + G. 


(b) 


In this table numbers refer to weight of air-dried material 
in gm and capital letters refer to the material; K = kiesel- 
guhr, S = maize starch, C = chalk, G = gap. Bracketed 
terms refer to layer increments, e.g. cake No. 4 was composed 
of an initial 300g starch, sith cakes formed by 100 g starch 
increments up to a total of 700g starch, followed by 100g 
kieselguhr layers up to a total of 300g kieselguhr. The test 
series thus included flow measurements on eight cakes, the 
final form being named cake No. 4. Cake No. 11 was composed 
of three layers only. The cake No. 20 included a gap between 
the chalk layers, whereas cake No. 22 had the gap support 
resting on the single chalk layer. 
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determining r,, Tp, and rye. Keeping the basket speed con- 


stant, another mass of slurry was added to form a new layer 
on the cake. By repeating the sequence, multi-layer cakes 
were made and their mean permeabilities determined. 
Layers were formed from two materials having a wide 
with a mean per- 


difference in permeability ; maize starch 


meability of about 3 x 10-7 g /sec.* and kieselguhr with K about 
4 x 10-% g /sec.? The radial position of the higher permeability 
layer was varied from being the outermost, innermost or middle 
layer between starch layers of different thicknesses. 

Table la specifies the different cakes. Cake No. 4 
was initially a 300g starch cake, then a layer of 
100 g starch was added, the new 400g cake being 
tested. Further cakes were made by 100 g additions 
up to 700g starch. This was followed by a layer of 
100 g of kieselguhr on the starch cake, the composite 
cake then being tested. This progressive addition of 
kieselguhr was followed up to a total of 300 g kiesel- 
guhr. After testing every final cake, samples were 
removed using an annular cutter, and the samples 
tested for filtration permeability in a 1} in, diameter 
filter cell [3]. Samples were cut at different positions 
of the cake and mean values of permeability are 
reported. The sample permeabilities were always 
within 2% of their mean. 

The overall mean hydroextraction permeability 
was calculated from : 

, 
o.g. log, r 
* Kem 


0 


* log, 


burn ‘o* 


C, log, (7 r) a 


where C, will be a constant for tests at constant yp, 
My Tor Tha» "he 

The mean filtration permeability was estimated 
from the following equation for the outflow time t, 
for volume V’ through a cake under a constant head 


H with inlet at a pressure slightly above atmospheric : 


uLV 


- 
AHK pp 8) 


A series of banded cakes was also formed and 
tested in a small filter cell which could be attached 
to the hydroextractor [3]. The layers were always 
formed by filtration under a suction head of 73 cm 
Hg and were tested for hydroextraction and filtration 
permeabilities. Table 1b specifies these cakes. In 
cakes Nos. 12, 13 and 14 a Perspex diaphragm was 
introduced to form a gap between the existing layer 
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and a new layer to be added. The diaphragm consisted 
of two dises, with perforations similar to the cell base, 
separated by three thin short rods. Two filter papers 
were introduced on the outer sides of the discs. With 
this diaphragm the effect could be studied of having 
a 1-2cm wide gap separating the cake into two 
layers. 

The materials used in the small cell experiments 
were chalk of approximate permeability 4-7 « 10-7 


9 


g/sec.? and kieselguhr with a permeability about ten 


times that of the chalk. During hydroextraction 








cake No symbol 
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Fig. 5a. Representation of equation (7) for homogeneous and 
banded cakes. 


tests on this cell, the cake surface at r=r, was 


protected by glass wool. }3] 


Discussion Or RESULTS 


The linear relation found (5) between the drainage 
time, t, from equation (5) and either (7,2 — r,*) or 
the mass of cake W,, still held for the multilayer cakes 
in the 9 in. diameter hydroextractor. Lines A and 
B in Fig. 5 represent total.starch or total kieselguhr 
cakes respectively. During the formation of banded 
cakes, the first layers of one material always followed 
the appropriate line, illustrating the reproducibility 
of the results. 
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In cake No. 4 (Table la and Fig. 5a), line .1 was 
reproduced up to the 700 g starch used ; then the part 
representing the kieselguhr layers was linear but not 
parallel to B. 

The test on the initial 200g kieselguhr layer in 
cake No. 3 is represented by point a on line B in 
Fig. 5a. When a thin layer of starch was added, 
point ¢ was obtained, followed by the other points 
for this cake as further starch layers were added. The 
line ae was steeper than the succeeding line, the 
decrease in the overall perméability being more than 
that expected from a calculated mean for 100 g starch 
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Fig. 5b. 


addition to a 200 g kieselguhir cake. Point a is based 
on the value of (r,? 
cake alone ({2], Fig. 8, curve B). If the layer is covered 
by a layer of starch, its thickness is reduced and a 
lower value of (r,2 —r*) obtained for the same 
kieselguhr mass ({2], Fig. 3, curve C). If the reduced 
(r,2 — r2) value for the compressed kieselguhr is taken 
as abscissa, the point a in Fig. 5a is moved to point 
b. The adjusted point b lies on the full line which is 
nearly parallel to line .4. The same phenomenon was 
noticed in cakes Nos. 5, 6, 7 and 8; the effect on the 
permeability represented by the ordinate in Fig. 5 of 


r2) from r, for a kieselguhr 





the first layer of starch added to the kiesclguhr cake 

was more than that expected for the increment in 

thickness represented by the abscissa, due to the 

compression effect associated with it. In each case 

the point b represents the point a moved to a radius 

equivalent to a kieselguhr layer below starch [2]. 
The condition represented by : 


t= C(r? —r?2 + Y) (7) 


was valid for additions within a band of one material, 
but the value of C varied from those for homogeneous 


Table 2. The value of © in equation (7). 


Cake No. Layer Cc. 

2 Kieselguhr alone 4-48 
Kieselguhr on starch : 

6 100 K on 3508 1-08 
5 200 K on 3508 1-92 
7 300 K on 3508 2-28 
s 200 K on 200 S 1-28 
4 300 K on 700 S 2-08 
1 Starch alone 0-448 


Starch on kieselguhr : 


3 700 S on 200 K 0-44 
8 500 S on 200 K on 2008 0-36 
6 350 S on 100 K on 3505S 0-40 
5 350 5S on 200 K on 350 S 0-312 
7 350 S on 300 K on 3508 0-232 


cakes of the particular layer material (Table II). From 
equations (5) and (7) 


» C, log, 2 

Are (8) 
whence an increase in C indicates an appropriate 
decrease in K. Thus kieselguhr exhibited a lower 
permeability when it formed a layer between starch 
layers than when it was alone. The decrease in 
permeability was greater for thinner kieselguhr layers 
than for thicker layers. 

All lines for the second layer of starch drawn from 
points 6 tend to meet at a point on line A where 
r, = 8-5 cm (Fig. 5), i.e. the overall mean permeabili- 
ties of a cake r, = 10-8 and r, = 8-5 was the same 
whether it was formed completely of starch or con- 
tained a middle layer of kieselguhr, regardless of 
either the position or the thickness of the kieselguhr 
layer. This suggests the necessity of a more detailed 
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study of these results to find out whether the middle 
layers of higher permeability were running full of 
water. 
The overall mean permeability of all tested cakes 
was less than that calculated as K,, from : 
" 


Pose 
1 log, nl 


r; ls - 

K, - kK, . in kK, 
For a cake such as the final cake No. 4 the general 
equation (9) applies with n = 2. The dimensions 
(r, — 1,) of the outer starch layer (700 g) and (r, — r,) 
of the kieselguhr layer (300 g) were measured. The 
permeability A, of the starch layer was measured 
after the cake had been spun at a speed higher than 
the test speed. This value corresponded to that shown 
on curve B of Fig. 10 in [2] and included the cloth- 
caek interface resistance R. The double-layer cake 
with 300g kieselguhr added was tested after com- 
pression at high speed and gave a mean permeability 
K,,, = 5-70 x 10-7 (Table 3). A permeability A, for 


Table 3. 


, 
log, 


(9) 


Comparison of estimated and measured mean permea- 
bilities for banded cakes. 





Mean Permeability, (K x 107 g /sec.*) 


K,, calculated from 
equation (9) 


Experimental 
Kom 


5°36 
5-93 
4-90 
4-04 
5°27 
4-90 





| 
(= 
| 
| 
| 
| 





the kieselguhr layer was assumed to be the same as 
that for a layer of the same dimensions and radial 
position in a cake of kieselguhr alone. This value was 
derived from mean data given in Table 4 of [2] by 
the application of equation (9). Using K, and K, in 
equation (9), the calculated mean permeability 
K,, = 5-93 x 10-7 was obtained. The slight differ- 
ence between K.,, and K,, justifies the assumption 
that the outer layer permeability remains constant 
when covered by inner layers, when all cakes are 
compressed initially at a speed higher than test 
speeds. 





For a three layer cake such as cake No. 6 the general 
equation (9) applies with n = 3. The dimensions 
To Ty» Tg Ts Were measured. Permeability K, of the 
initial layer was measured and assumed constant 
thereafter. Permeabilities A, and K, were deduced 
from previous data [2] as described above. The 
calculated mean K,, and the measured permeability 
K_,, for the banded cakes are given in Table III. 

The differences between K,, and K.,, for cakes 3, 
5, 6, 7, 8 may indicate that either the banded cakes 
offered more resistance than expected from equation 
(9), or that the flow through the kieselguhr layer was 
discontinuous and hence the centrifugal head deve- 
loped across them was less than that to be expected 
from equation (3). 
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Figs. 6, 7, 8. Pressure distribut'on within banded cakes in the 
hydroextractor, with n’ = 1,050 r.p.m. and ry, = 8-0 cm. 





8 


Using the layer permeabilities K,, K, derived as 
for insertion in equation (9) and the measured per- 
meability K, of the first layer of a banded cake, the 
friction head P, from r =r, to r =r at any radial 
position in the cake could be found from equation (4). 
The values are represented in Figs. 6, 7 and 8. The 
P. line A gives the centrifugal head developed up to 
any radius within the liquid, when r,; = 8-0 cm. 
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Taking cake No. 5a (Fig. 6) as an example, the point 
a represents the head developed up to the back of 
the innermost starch layer; the vertical line ab 
represents the frictional head, P,, dissipated across 
this inner layer of starch. The line be represents the 
centrifugal head developed in the kieselguhr layer and 
is parallel to the P, line for the same radial positions, 
assuming that this layer was running full. As the 
pressure at the back of the cake at r = r, is assumed 
to be atmospheric, so fg will represent P, across the 
outer starch layer. The variation of P, given by ef 
must be parallel to line A, assuming that this starch 
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layer runs full, e being at the inner radius of the 
outermost layer. 

The pressure drop across the kieselguhr estimated 
from K, and the known flow rate q is ed. If all three 
layers were running full the friction heads should sum 
to ab + ed + fg = hg. The head hg is not accounted 
for, a considerable discrepancy cd being found in 
each case. The explanation is presumed to be that 
the cake is not running full at all radii. Cakes of one 
material conform to equations based on the con- 
tinuous development of P. in line A at all positions 
(Fig. 1, 2, 3), the outermost layer resistance is enhanced 





by the effect of the cloth interface, and the kieselguhr 
is ten times as permeable as starch. The kieselguhr 
layer was probably not running full and was thus 
not developing P, comparable to line A, while the 
void space in both starch layers was completely 
filled with liquid. The development of a head bd 
within the kieselguhr layer is discussed later together 
with the results in Fig. 14. The same phenomenon 
was noticed in all cakes where a kieselguhr layer 
existed between layers of starch (Figs. 6, 7 and 8). 
The discrepancy cd was larger with thicker layers of 
kieselguhr (Fig. 7). 
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For cake No. 4, where kieselguhr formed the 
innermost layer, the calculated friction head added 
together to approximately the centrifugal head 
developed, giving very small discrepancy cd (Fig. 8). 
When the kieselguhr formed the outer layer (cake 
No. 3), the pressure at the outer surface of the starch 
was atmospheric, indicating that the kieselguhr layer 
was not running full of liquid (Fig. 8). 

Samples were cut from these cakes and tested for 
filtration, giving the results in Table 4. It was not 
possible to cut samples from cakes made from kiesel- 
guhr alone. It was assumed that the mean filtration 
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Table 4. Filtration tests on samples of banded cakes from hydro- 





Kp x 107 
for kieselguhr 








permeability for starch of 3-75 x 10-7 g/sec.? found 
from cake No. 1 was the mean permeability in all 
the other starch layers in succeeding cakes. The 
permeability of the kieselguhr layer given in the last 
column of Table 4 was calculated from the relation 


a "n (10) 


, 








permeability K, = 46-0 x 10-7 for inner layers of 
kieselguhr in the hydroextractor. This suggests that 
in filtration all layers of banded cakes run full, except 
the bottom layer, if this is of highly permeable 
material. In the case of cake No. 3, the calculated K , 
for kieselguhr was negative, showing that the per- 
meability assigned to the starch layer was low. When 
the mean K, of inner layers of starch (K, = 4-6 x 10-7) 
was assumed to represent K, for the starch layer, Kp 
for kieselguhr was found to be infinite, the kieselguhr 
layer in cake No. 3 apparently offering no resistance 
to flow. This is consistent with the filtration head 
being dissipated across the starch layer, the liquid 
only running through some of the pores of the kiesel- 
guhr under the gravitational field. This supports the 
earlier assumption [3] that the highly permeable cloth 
runs empty in all the hydroextraction and filtration 
tests. The result also confirmed the existence of a 
resistance at the cloth-cake interface. The mean 
filtration permeability was 3-75 x 10-7 whenever the 
starch was adjacent to cloth, for starch layers not 
in contact with cloth the mean was 4-6 x 10-7 g/sec.” 
The difference was due to the inclusion of the cloth- 
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Fig. 9. Relation between rate of flow and speed of revolution for banded cakes. Line N is for the 5 g chalk cake data in A. 


115 





M. M. Harun and J. ANDERSON Storrow : Hydroextraction V : Radial distribution of pressure in homogeneous cakes 


starch interface resistance R in the nominal resistance 
of the starch layer adjacent to cloth. 


Tests IN THE SMALL Spinninc CELL 
When cakes in Table 1b were tested for the charac- 
teristic g, n data [3] it was found that banded cakes 
still followed the relation derived for a homogeneous 
cake : 


(11) 
The logarithmic g, n plot was a straight line of slope 


nearly 2-0 for all banded cakes (Fig. 9) and the overall 
mean permeability was independent of speed. The 























Litiliyyy 


° $00 
TIME ,t seconds 





1000 


Fig. 10. The application of equation (12) in testing representa- 
tive banded and homogeneous cakes. 


plot of log H against ¢, for filtration under a falling 
positive head with discharge at atmospheric pressure 
was linear (Fig. 10), satisfying again the filtration 
equation (12) derived for homogeneous cakes, for the 
time t, required for the head to fall from H, to H, 
in a stand-pipe of cross sectional area a 
H, K,.A.t, 
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Figs. 11, 12. Pressure distribution within banded cakes in the 
small spinning cell, with n’ = 1,500 r.p.m,. and rz, = 0. 
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From hydroextraction tests on the chalk and kiesel- 
guhr used in these experiments, the layer permeability 
and cloth-cake interface resistance were calculated. 
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The friction head across each of the layers in a 
banded cake was found and plotted in Figs. 11, 12. 
From Fig. 11, the mean discrepancy cd forms 3.4% 
of the total P., which could be explained on the 
assumption that the highly permeable middle layer 
was not running full. When kieselguhr formed the 
outer layer (Fig. 12) the pressure at the back of the 
chalk was atmospheric giving another example of a 
layer of kieselguhr running empty. The calculated 
friction loss across such layers de is much less than 
the expected centrifugal force ce, which could develop 
if the outer layer was running full of liquid. The 
discrepancy is again greater for thicker layers of 
kieselguhr. 

The filtration tests for equation (12) were conducted 
here by subjecting the cake to a falling positive head 
around 120 cm water, without applying any suction 
under the cake. The time t, for the head to fall 
from H, to H, was found and the mean overall 
filtration permeabilities calculated from equation (12) 


Table 5. Filtration tests on banded-cakes in small spinning cell. 





Kp x 107 
Jor kieselguhr layer 


5 g Chalk 
11 
12 
13 
14 
15 
16 
17 
18 
19 








are shown in Table 5. The permeability of the 
kieselguhr layer was estimated from K, = 4-9 x 107 
for chalk alone and are given in the last column in 
Table 5. Apart from the expected differences due to 
compression stresses of successive layers, the results 
indicate that the highly porous layer was still running 
full, even when it formed the layer of the cake dis- 
charging to atmosphere. 

In the cakes Nos. 20, 21, 22, the light Perspex gap 
support did not cause any appreciable compression 





when resting on the 5 g chalk cake, when tested for 
q, n data (Fig. 13). When the gap was used to separate 
the chalk into two layers as in cakes Nos. 20, 21, line 
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Fig. 13. Relation between flow rate and speed of revolution 
for cakes separated by a gap. 





B in Fig. 18 was obtained. After the test was com- 
pleted, the liquid feed was stopped and the cake was 
allowed to run dry for a minute, followed by a rapid 
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Fig. 14. Pressure distribution within cakes separated by a gap. 


117 





M. M. Haruni and J. ANDERSON Storrow : Hydroextraction V : Radial distribution of pressure in homogeneous cakes 


test of the g, n condition. With cake No. 21 no 
difference could be observed between the first rapid 
test and others taken later, all agreeing with line B. 
However, in the case of cake No. 20, the first point 
found quickly after the feed was resumed fell on line 
C, whereas the other points determined afterwards 
returned to line B. 

When the friction drops across the layers were 
calculated and plotted in Fig. 14, the gap resistance 
de being zero, it was found that cd formed 17-5% of 
the total expected centrifugal head hg. The centrifugal 
head fd was extrapolated to a point j at the same 
pressure level as 6. An empty compartment of 
constant pressure probably extends behind the inner 
layer of cake from the radial position b to position j. 
The section of the gap from j to d was filled with a 
layer of liquid maintained in front of the outermost 
chalk layer. The friction head dissipated in this 
outer layer was the sum of the centrifugal head 
developed from j to d in front of the cake, from d to 
f within the cake, and the pressure differential between 
the empty compartment and the atmosphere, repre- 
sented by the ordinate to b. 

When the cake was allowed to run dry, the pressure 
in the gap decreased to atmospheric, as the outer 
layer emptied, then began to increase when the flow 
was resumed. 

In the case of cake No. 21, the inner layer would 
allow a flow of 1-58 cc/sec. when running full and 
discharging to atmosphere at the radial position }, 
whereas the outer layer would allow 0-22 cc ‘sec. when 
a liquid layer extended from position ) to r = ry 
at g. Due to this great difference, the pressure in the 
gap would be increased on the resumption of flow 
and readings showing lower flow rates than line B in 
Fig. 13 could not be detected. The corresponding 
values for cake No. 20 are 0-66 and 0-31 cc/sec. 
respectively, giving an opportunity for a rapid test 
to detect lower flow rates, as on line C, Fig. 13, when 
the pressure in the gap was still developing. 

Filtration results conforming to equation (12) 
provided filtration permeabilities given in Table 6. 
The pressure of the gap support and its filter papers 
in cake No. 22 had an insignificant effect on permea- 
bility. The 6% reduction in permeability with two 
chalk layers is probably associated with an interface 
resistance between chalk and filter paper, mainly at 
the back of the upper chalk layer. The great effect 


Table 6. Filtration tests in small spinning cell. 





Cakes with layers separated by a gap 


Cake No. | Kp x 107 
— = on | - —EE —— 
5 g Chalk 49 
20 4:57 
21 4-66 
22 | 4-83 








noted in hydroextraction attributed to the gap 
running empty did not occur in filtration when the 
filtrate was discharged at atmospheric pressure. 


CONCLUSIONS 


The study of hydroextraction and filtration per- 
meabilities of banded cakes has demonstrated that 

(i) a backing layer of high permeability material 
in a hydroextractor does not run full of liquid and 
the atmospheric pressure is exerted on the outer face 
of the inner layer of less permeable material. This 
conclusion based on 10: 1 ratio of permeabilities is 
considered to confirm the assumption that a cloth 
support in the basket does not run full, the corres- 
ponding permeability ratio being 50 : 1. 

(ii) the given interpretation of the results is in 
accordance with the previous conclusion that there 
is a small effect of compression on cake permeability 
with the present cake materials, and a significant 
resistance to flow at the cloth-cake interface, 

(iii) a highly permeable mid-layer within a cake 
does not run full of liquid, a constant pressure region 
existing within this layer, 

(iv) the conditions here demonstrated to cause a 
flow discontinuity within a cake are absent with 
cakes of one material, and the assumption of con- 
tinuous centrifugal head development within homo- 
geneous cakes, fundamental to the equations and 
proposed mechanism of hydroextraction, is supported, 

(v) with usual operating conditions, it is improbable 
that a suction head develops in any region within the 
cake to form a discontinuity in the assumed flow 
pattern. 
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NOMENCLATURE 


= cross-sectional area of filter bed 

= cross-sectional area of tube in filtration tests 

== gradient in equation (7) 

= gravitational acceleration 

= hydrostatic head in filtration tests 
H,, Hg limits of 11. 

= permeability ; K,, Kg,...K, for different 
layers ; K p for filtration ; K, for hydro- 
extraction ; K pm, K-» appropriate experi- 
mental means ; K,, calculated mean from 
K,, Kg, ... Ky 

= length of filter bed in filtration tests 

= speed of revolution 


= centrifugal pressure 

= friction head 

= drainage rate through cake 

= cake-cloth interface resistance 


cm? 

em? 

sec. jem? 
em. /sec.? 
cm 


g /sec.* 
em 

rev. /sec. 
rev. /min. 
cm water 
em water 
ce /sec 
sec.? /g 





= radius from basket axis (cm); rg for cake 
outer surface ; r, for cake inner surface ; 
r,, for liquid inner surface ; 7), Tg. - . - Tn 
for interfaces between cake layers ; rz), 
Tr limits of r; in drainage time equation 

= time ty, time in filtration equation 

= volume flowing in time {tp in filtration 
equation 
mass of dry solid in cake 

= height of centrifuge basket 

viscosity of liquid 
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Summary—The conditions of validity of the simple forms of the Gibbs-Duhem equation 
most commonly used in practice are discussed. A simple derivation is given of the exact relation- 
ships valid at p constant (resp. T constant) in the case of a binary vapour-liquid system at 
equilibrium and connecting together either (a) the fugacities of the components, or (b) the activity 
coefficients of the components, or (c) the compositions of the vapour and liquid phase with T 
(resp. p). The application of the general relationships is discussed, especially for the case of a 
constant pressure, important in distillation problems. 


Résumé—Les auteurs discutent les conditions de validité des formes approchées de I’équation 
de Gibbs-Duhem employées dans la pratique. Sont également discutées les équations exactes 
valables soit & p constant, soit 4 T constant, dans le cas d'un systéme binaire liquide-vapeur en 


équilibre et reliant entre elles soit: (a) les fugacités des constituants ; (b) les coefficients 
d’activité des constituants ; (c) les titres des phases gazeuse et liquide en fonction de T ou de p, 


respectivement. 


The Duhem (or Gibbs-Duhem) equation is 
commonly applied to binary, liquid-vapour 
systems at equilibrium and is widely used in 
the calculation of the composition of the coexisting 
phases of binary systems important in distillation 
problems. Two of the most often encountered 
forms of this equation are 


din f, _{ ding, 
(4) " (sean a).. (1) 
and 
dIn y, a d In y. 
( ding ). ' (; In (1 a) (2) 


The derivation of these equations, as given for 
instance by Lewis and Ranpa.t [5], restricts their 
validity to conditions of constant pressure and 
temperature. In the case of a two-phase system 
these equations are therefore no longer equilibrium 
relations, since the phase rule shows that in such 
a system one cannot have at variable composition 
both temperature and pressure kept constant and 
have the equilibrium between the two phases still 
being maintained. In vapour-liquid calculations, 
therefore, one of the two conditions (p = const. 
or T' = const.) has to be dropped ; eqs. (1) and 





(2) are then used as strictly isothermal ones, 
whereas eq. (2) has also been used by Beatty 
and CALINGAERT [1] and others at constant p 
and variable T. Both ways of using eqs. (1) and 
(2) are only approximations and it is important 
in a given case to discover if these approximations 
are admissible or not. This is especially important 
when eq. (2) is used as an isobaric relation, because 
then the approximation is often not very good. 
The object of the present note is to give a review 
of the general and exact relationships valid at T 
constant and p variable on the one hand, p con- 
stant and T variable on the other, the relations 
being derived in a simple but rigorous fashion 
which shows clearly the conditions of validity 
of the general forms of the Gibbs-Duhem equation 
as well as that of the more simple relations com- 
monly used in practice. 

The thermodynamic potential F of a homogen- 
eous system is a homogeneous function of the 
first degree in the mole numbers of the compo- 
nents of the system. From Euler’s theorem on 
homogeneous functions follows (for a binary 
system) 

F = 1 py + Mg Me 
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and therefore 
(3) 


On the other hand we have the familiar relation, 
valid for any open phase : 


dF =n, du, +n, dp, + pw, dn, + pw. dn, 


dF = Vdp — SdT + p, dn, + py dn, (4) 
Subtracting (4) from (3) we get 
Vdp — SdT = n, dp, + n, dp, 


This is the most general form of the Gibbs equation 
valid at variable T and p. It may be written as 
well for the liquid phase as for the vapour phase 
of a binary two-phase system at equilibrium. 

If we divide by the sum of the mole numbers 
n, +n, and bear in mind that the sum of the 
mole fractions is 1 we get, for instance for the 
liquid phase 

V~,dp —S,dT = adp, + (1 — 2) du,t (5) 
At constant temperature and pressure this 
equation reduces at once to eq. (1), if we take 
into account the relation between » and the 
fugacities 


(6) 


(a) Relationships valid at T constant and p variable 


At constant temperature eq. (5) takes the form 
(if we bear in mind the relation (6)) 


at T constant. 


du = RTdlnf 


ae 
RT 


If we regard z, the mole fraction of component 1, 
as the main variable we may also write 


dinf, dinf, 
ding din(1l —2) 


adinf, +(1—a)dInf,=—"dp at T constant (7) 


V., dp 


RT ar at T constant (8) 


+ 


or 

Vi, dp 
ding din(li—a2) | RT dr 
the latter equation applying when the gas phase 


dinp,___ dinp, 


at T constant (9) 


+ We do not give a superscript to the variables p,, ys 
p and T since at equilibrium they have the same value in 
both phases. 





is ideal. Eqs. (7) and (8) are exact forms of the 
Duhem equation which remain strictly valid if 
we apply them to equilibrium calculations in a 
two-phase system maintained at constant tem- 
perature. Eq. (8) differs from the simple relation 


;, Vm dp dp 
(1) by the term RT ar’ where y 
the curve of the total vapour pressure vs. the 
composition of the liquid phase (at constant 7’), 
and V,, is the volume of 1 mole of the liquid 
phase, which may easily be calculated from the 
density of the solution. It will be possible there- 
fore in a given case to evaluate at least approxi- 
V, dp 
RT as and thus 
the degree of accuracy given by the simple 
relation (1). In general, this latter can be used 
without significant error if, as is usually the 
case, the molal volume of the liquid phase is 
negligible in comparison to that of the vapour 
zdlnp, 


is the slope of 


mately the correction term 


phase (this because of the fact that and 


— 2x) d In p, 
(1 lh MP2 are larger or at least of the same 
( 


iln p 


; c —- 
order of magnitude as —_“ , which is equal to 


fo) - in the case of an ideal gas phase). 

Equations (8) and (9) give a relation between 
the partial pressures (resp. the fugacities) of the 
two components of the mixture. Often, however, 
one wants a direct relationship between the total 
pressure p and the compositions z and y of the 
liquid phase and the vapour phase (that is an 
equation relating the dew-point line to the bubble- 
point line). The general form of this p-a-y 
relationship may be derived as follows. 

The fugacity f, of the component 1 of the vapour 
phase is a single-valued function of T, p and y 
(or Iny). We may therefore calculate the total 
exact differential of In f, as follows 


OA) or +f) ans 


dInf, 
+ (3 = A) diny (10) 


dinf, = 


t See also footnote, p. 122, column 2. 
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An analogous relation holds for f, : 


ante = (Sar) s (Sp')ees 


dinf, ) 
anal idin(l — 11 
+ (525 a ame y) (11) 
We have also the familiar relations (see for 
instance [2d], [3b] or [4b)). 


C34) = a (2) 
?24) =a (18 


where }, and 0, mean the partial molal volumes 

of component 1 and 2, respectively. 
Furthermore, we may write on the basis of 

(5) and (6) applied to a single phase at constant 


T and p 
(4)- (eae ode . (14) 


At equilibrium the fugacities are the same in both 
phases. We may therefore replace in eq. (7) 
dinf, and dinf, by their values as given by 
eqs. (10) and (11) which readily yields (if we bear 
in mind eqs. (14), (12) and (13) and the fact that 
the first term of the right-hand side of eqs. (10) 
and (11) is zero at 7 constant) : 


a Cr) at T constant (15 
RT y(l—y)AV \diny “yA (15) 
where AV is equal to 2 0, + (1 —2) i, — V,; 
it is the increase in volume when 1 mole of the 
liquid phase is evaporated into such an amount 
of the equilibrium vapour phase that the com- 
position of the latter is not substantially changed. 
The partial derivative (F) T refers to the 
diny T.p 
vapour phase and means the variation of fugacity 
of component 1 with composition at constant T 
and p. Eq. (15) is the exact p-a-y relationship 


t As has been already pointed out the partial derivative 


(34) does not correspond to changes along the 
diny T.p 

equilibrium line, when the vapour phase is at equilibrium 
with the liquid phase. It still has, of course, even in 
this case a perfectly determined value. 





valid at T constant. If the vapour phase is a 
mixture of ideal gases we have 


4) - a BE 
diny T.p 


If then the volume of the liquid phase is negligibly 
small, eq. (15) reduces to 


dp me ee dy at T’ constant (16) 
Pp y(l—y) 
This equation, as well as 
(22%) “ (pa Ps (17) 
ding T din(1 — 2) T 


are two of the forms of Duhem’s equation, which 
are most commonly used in practice. They are 
isothermal relations valid only if : (1) the vapour 
phase may be considered as a system of ideal 
gases, (2) the molal volume of the liquid phase 
is negligible in comparison to that of the vapour 
phase.t If these conditions are not fulfilled the 
exact relationships (15) or (8) should be used. 
The use of the exact eq. (15) is quite similar to 
that of eq. (16) (see for instance [2c]). The values 


of 0, d, and (° =) have to be taken 
T.p 


from an equation of state of the real gas-mixture. 
A considerable simplification occurs if we admit as 
a first approximation the validity of the Lewis- 
Randall rule 


j= uh f=(l- y) fe 
(ft. fe = fugacities of the pure components) 


+ It should be noted that this latter condition does not 
hold quite in the same way for eq. (16) as for eq. (17). In 
the case of eq. (16) the relative volume of the liquid phase 
corresponds exactly to the error involved when using this 
equation, that is, if the molal volume of the liquid phase 
is for instance 0-1 % of that of the vapour phase, both sides 
of eq. (16) will also differ by 0-1%. Such a simple relation 
does not hold in the case of eq. (17). In general if the 
relative volume of the liquid phase is for instance 0-1% 
both sides of eq. (17) will also differ by less than 0-1% or 
0-2%. But if the vapour pressure of one of the pure 
components of the mixture is much larger (for instance 
four times larger) than that of the other component a sub- 
stantially larger error may occur at one of the ends of the 
diagram, where z or |-z is small. 
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which holds in general with a fair approximation 
at pseudoreduced pressures less than 0-8 [4c]. 
( Inf, 

diny 
are equal to the molal volumes of the pure com- 
ponents in the gaseous state at the temperature 
and pressure of the mixture. They may be 
evaluated from an equation of state of the pure 
components (or from a generalized (reduced) 
compressibility diagram).f 


is now equal to 1 and 0, and i, 
T.p 


(b) Relationships valid at p constant and T variable 
As most vaporization operations are conducted 
at constant pressure and variable temperature 
this case is of considerable practical importance. 

If we write the analogues of eqs. (10) and (11) 
for the liquid phase we have 


ands = Cape) ae + ar) + 


4 pA) ar ding (18) 


ses Oa CHB) ar 


Jaina —a) (19) 
T.p 


the first term of the right-hand sides of these 
equations being zero at p constant. Furthermore 
we have the relations (see for instance [8c]) 


ei), 


O),. 


where H; and H;, are the molal enthalpies of the 
pure components in the ideal gas state. 


_H;-; 


RT? ” 


+ It may happen that one (or both) of the pure gaseous 
components would be liquefied if it existed alone at the 
temperature and pressure of the mixture. The equation 
of state (or the reduced compressibility diagram) has then 
to be extrapolated into the two-phase region. This pro- 
cedure will not be very accurate if extended extrapolation 
is necessary. In general, however, only a rough knowledge 
of the departure of the vapour phase from ideal behaviour 
is needed in the use of Duhem’s equation. 





Combining (20) and (21) we may also write 


(nh) 4 a —2) (4) AH 


oT MT ) 1-2 RT? 


where AH means the term 2H; +(1—2) H, —zH, 
— (1 — 2) H, and is therefore the molal ideal heat 
of vaporization (i.e. the heat required to vaporize 
1 mole of the solution, when the pressure of the 
obtained vapour is sufficiently low for the ideal 
gas laws to be valid). 


(22) 


Rearranging eqs. (18) and (19) and combining 
them with eq. (22) we readily get, if we bear in 
mind that ( : 7A) and (n* 5) are 

IInz/r» din (1 —2)/7> 
equal because of Duhem’s equation when applied 
to a single phase at T' and p constant 


dinf,  dinf, AH dT 


= —— t 3 
ding din(i—a)* CemeaeR 





RT? dz 
Eq. (23) is an exact form of Duhem’s equation 
valid at constant total pressure and variable 


temperature. - is the slope of the curve T vs. 


x at p constant (boiling point line). AH may be 
evaluated from the heat of mixing (corresponding 
to the composition z of the solution at the con- 
sidered point of the diagram) and the ideal heats 
of vaporization of the pure components. In most 
cases the heat of vaporization is large in com- 
parison to the heat of mixing, so that this latter 
needs often to be known only roughly or may even 
be neglected. 

It should be noted that at constant pressure 
the simple eq. (1) cannot be used even as an 
approximation. This can be seen in the following 
way. 

If we assume the gas-phase to be ideal we have 


fe =(1—y)P 


AH dT 
RT? de eq. (23) may 


Sf, = yp 
If we then neglect the term 


be put into the form 


o-s 
y (1 —y) 

which relation obviously does not hold even 

approximately. 
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Eq. (23) relates the curves of the fugacities — or 
in the case of an ideal gas phase, of the partial 
vapour pressures — of the two components with the 
composition of the liquid mixture, when the total 
pressure is kept constant. In isobaric problems, 
however, equations relating directly the boiling 
point line to the dew-point line are often more 
useful. The exact T7-a-y relationship may be 
simply obtained as follows. At p constant the 
second term of the right-hand side of eqs. (10) and 
(11) is zero. If we replace in eq. (23) d In f, and 
din f, by the values given by eqs. (10) and (11), 
noting that (° 2) and ( a Inf, 

Iiny/ 7,» dIn (1 —y)/ 7,» 
are equal, we get 


AH a) dinf, 
— Tf —(1— iT = 
ler (Sr), Cor al 
= (7) ee dy atpeonstant (24) 
diny/r,py(l — y) 


The terms C4") ; ems) 


oT dT 
dinf, 
( d In y 
are the corrections for the imperfection of the 
vapour; the first two are equal to zero and the 
third to unity when the vapour phase is a mixture 


of ideal gases. In this case eq. (24) reduces to 
e#—y 


oa = y(1—») AH dy atpceonstant (25) 
which is at p constant the analogue of eq. (16). 
Eqs. (25) or (24) may be used to check the 
consistency of experimentally determined 
vapour-liquid boundary curves in a way quite 
similar to that ¢ employed when applying eq. 
(16). The most accurate method is to integrate 
to determine T (the righthand side of eq. (25), 
and in the case of eq. (24) both sides of the 
equation, have to be integrated graphically).tt A 
further possibility is to calculate the value of 


and 
p.i-y 


refer to the vapour phase: they 
Tp 


“om eq. (25) or (24) and compare it with the 


slope of the dew-point line. Other possible 
applications of eqs. (24) and (25) are the calculation 
- t The use of eq. (16) is given for instance by Donor [2c]. 

tt The way of evaluating AH has already been dis- 
cussed. 





of the boiling-point line, if only the dew-point 
line is known, or the evaluation of the heat of 
mixing (and thus of the enthalpy of the liquid 
mixture) from the T-2-y diagram and the heats 
of vaporization of the pure components. 

When the departure of the gas phase from 
ideal behaviour cannot be neglected, eq. (24) has 
to be used. The correction terms for the imper- 
fection of the vapour have to be evaluated from 
an equation of state of the real gas mixture. A 
considerable simplification arises if we assume 
the validity of the Lewis-Randall rule. (; 2) 

, IIny/7,» 
is then equal to 1 and the two other correction 


terms are given by 


dinf, H; —H;° dInf, H,—H;, 
Pf), -MatiaaA) - Hat 
oT / py oT} p,1-y RT? 
where H; and Hy are the molal enthalpies of 
the pure components in the ideal gas-state whereas 
H,* and H,° are the molal enthalpies of the 
pure gaseous components at the temperature and 
pressure of the mixture. These correction terms 
may be calculated in the usual way from the 
equations of state of the pure components (see 
for instance [2a], [3a], [4a]), or from the generalized 
(reduced) enthalpy-diagrams (see for instance 


[2b]).f 


(c) General forms of Duhem’s equation in terms of 
activity coefficients 

If we choose as the standard states the pure 

liquid components at the pressure and tem- 

perature of the solution, the activity coefficients 

are given by 


, =f ia 

" wh (1—2) f.* 
where f,* and f,* are the fugacities of the pure 
liquids at the same pressure and temperature as 
the mixture. If we introduce these values into 
equation (8) we readily obtain 
adinffdp diny, 


ding © dp 
_(l —«) dln fy dp + 
dp dx 


t See also footnote, p. 123. 


Ye 


din y, 


dz din(i—z) 
V,, ap 


RT dz at T constant 
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"aad — — ‘tan ye 1 dp 
+|Pa aV, (1 0) ¥s'| gpa 


at T constant (26) 


where V,° and V,.° mean the molal volumes of 
the pure liquids at the pressure and the tem- 
perature of the mixture, and the term V,, —aV ,° — 
(1—2)V,° is therefore the volume change accom- 
panying the mixing (per mole of mixture). We 
will denote the second term on the right-hand 
side of (26) by the letter A. 
Similarly, we obtain from eq. (23) 


diny, ,adinffdT _ diny,  _ 
ding dT de din(l —2) 

1—az)dinfg dT . 
= © 1K. as + leu; + (1 — 2) WH, — 
dT 


at p constant 
daz \ 


' — 1 
H, —(1 —2) Hy} - 
r 1 ( 2) | RT? 
or 
diny, ____diny, 
ding din(1l —2) 


1 dT 


RT? da 
(27) 


- len, +(1—2) Hj —2 H;*—(1—2) us; 


at p constant 


where the term zi, + (1 —a2) H, —2H,* — 
(1 —a) H,* is the integral heat of mixing (per 
mole of mixture).| We will denote the second 


+ Strictly speaking this heat of mixing (as well as the 
volume effect of the mixing) should be measured at the 
pressure and the temperature corresponding to the com- 
position z of the considered point of the curves y, and y, 
vs. 2. The available experimental values will however in 
general correspond to another pressure than that of the 
considered mixture (this also because of the fact that the 
standard state as defined above may be a hypothetical 
one). This will however involve no appreciable error in 
the evaluation of the correction terms A and B. As these 
last values will need most often to be known only approxi- 
mately it will often even be sufficient to know the heat of 
mixing and volume change at another temperature than 
that of the mixture (for instance at room temperature). 





term on the right-hand side of (27) by the letter B. 

(26) and (27) are the exact equations relating 
together at 7’ and p constant, respectively, the 
curves of y, and y, vs. z. The approximate 
equation (2) may therefore be used (a) at T con- 
stant when the correction term A is negligible, 
(b) at p constant when the term B is negligibly 


In y, 


small in comparison to poh . The first condition 
ding 


is mostly fulfilled (except perhaps at high tem- 
dp 
dx 
may be large) and eq. (2) will in general hold at T 
constant with a high degree of accuracy. It is 
not so, however, at p constant. In this case before 
using the simple relation (2) the term B (which 
depends on the steepness of the boiling point 
curve and the heat of mixing) should always be 
evaluated at least approximately and compared 
to d Iny, 

dina 
the simple relation (2) might hold fairly well in 
certain parts of the diagram and still be in 
serious error at other values of z. 

Let us take as example the system acetone- 
carbon disulfide studied by various authors [6], 
[7], [8] at t = 35-2° and p = 760 mm (Beatty 
and CaLINGAERT [i] have applied eq. (2) to this 
system, to check the values obtained at T constant 
as well as those measured at p constant). At 


peratures and pressures where — and the volume 


It should be noted that at p constant 


Ip . 
z= 0-5, 5 is about 0-13 atm., and the volume 
c 


changes by about 1-4% when the components 
are mixed [7], which gives for A a value of 
din y, 
ding 
at T constant both sides of eq. (2) differ only by 
0-0007%. Again at z = 0-5 the heat of mixing is 


5 x 10°, being ca. 0-7, it follows that 


330 cal/mole [7], td is about 3-5° and B is then 
C 


equal to 0-005. At p constant both sides of eq. (2) 
differ therefore by 0-7% when 2 =0-5. The 
same calculation for z = 0-9 and 2 = 0-1 reveals 
that the degree of accuracy given by eq. (2) is 
now rather poor since B is not less than 7% and 

din y, 


18% of _.— 


dine’ respectively. 
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NOTATION 
= fugacity 
= thermodynamic potential (= H TS) 
= enthalpy 


partial molal enthalpy 

ideal heat of vaporization (per mole of mixture) 
number of moles 

total pressure 


= partial pressure 


ideal gas constant 

total entropy of a phase 

molal entropy (entropy of a mole of the mixture) 
thermodynamic temperature 


= total volume of a phase 


molal volume (volume of a mole of the mixture) 
partial molal volume 


= mole fraction of component 1 in the liquid phase 
= mole fraction of component 1 in the vapour phase 
= chemical potential (partial molal thermodynamic 


potential) 





y = activity coefficient 

Subscripts : 1 and 2 refer to component 1 and 2, respec- 
tively. 

Superscripts : ‘and “’ refer to the liquid and vapour phase, 
respectively ; ° refers to the standard state of low 
pressure, where gases are ideal; * refers to a pure 
component. 
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Summary— The accepted methods for the design of condensers to handle gas-vapour mixtures 
are reviewed and a single verticle tube condenser is described from which data are obtained for 
the systems chloroform and air, and steam and air. 

The method of CoLsurn and Houcen [4] is applied to a run on each of the systems and the 
calculated area compared with the expected area. The approximate method of Co_aurRn [3] is 
also applied to the same runs and the area obtained again compared with the expected area. 

The area required to cool a saturated mixture of chloroform and air from 120-2°F to 94-0°F 
calculated by the Cuo_Lpurn and HoucGen [4] method is found to be in good agreement with the 
expected area, being within 2%, of it. The method applied to the cooling of a saturated mixture 
of steam and air from 168-3°F to 129-9°F gives an area which also agrees well, being within 5% 
of that expected. 

The modification of the CocpurRN and HovucGen [4] method suggested by Sarru [10] is found 
unnecessary in the cases studied. 

The approximate design procedure proposed by CoLBuRN [3], applied to the same conditions, 
is found to give an area approximately 40% greater than expected for both systems. 


Résumé —L ‘auteur passe en revue les méthodes classiques de calcul pour les condenseurs traitant 
des mélanges de gaz et vapeurs. II décrit un condenseur a tube vertical unique utilisé pour 
l'étude des systémes air-chloroforme et air-vapeur d'eau: Les résultats expérimentaux sont 
comparés avec ceux prédits par les méthodes classiques. 

Dans le cas de CHC], —air, (condensation entre 120-2°F et 94-0°F), la méthode de CoLpurn — 
Hovcen [4] concorde & 20% prés. Dans le cas H,O—air (condensation entre 168-3°F et 129-9°F) 
la concordance atteint environ 5%. 

La modification apportée par Smiru [10] & la méthode précédente n'est pas nécessaire dans 
les deux cas étudiés. 

La méthode de calcul approximatif proposée par CoLsurn [3] prévoit les surfaces en excés 


de 40%, dans les deux cas étudiés. 


INTRODUCTION unit time per unit cross-sectional area through 
the resistances of the condensate layer, the metal 
wall, scale and cooling water film is made equal 
to the heat flow through the gas film. The heat 
flow through the gas film is made up of sensible 
heat lost by the mixture and latent heat trans- 
ferred as the vapour diffuses through the gas film 
and condenses in the condensate film. The total 
heat flow is represented by an overall coefficient 
U, multiplied by the temperature drop from the 
mixture to the water. 
These conditions are represented by : 


When a gas-vapour mixture is exposed to a 
surface at a lower temperature than the dew 
point of the mixture, condensation of the vapour 
occurs. A gas film is formed at the surface, 
which contains a higher gas concentration than 
the main stream. 

CoLBuRN and HoucGen [4] have shown that the 
diffusion of the vapour molecules through the gas 
film plays an important part in the transfer of 
heat from the gas-vapour mixture. They pre- 
sented a design procedure for estimating the 
surface area required for a cooler-condenser, to h, (t, —t,) + K M,A(p, — p,) 
cool a saturated gas-vapour mixture, along the = h, (t, —t,) = UAt (1) 
saturation line. sais 

In their treatment, at any particular point in| or in words, the heat flow to the condensate 
the condenser, the quantity of heat flowing per! surface is equal to the heat flow from the 
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(‘: — 'e) = 0-2 (3) 


condensate surface and these are each equal to 
U At. 

At any value of t, in the condenser, all the 
variables in equation (1) are known or can be 
calculated, except t,, p, and U. The value of 
p, follows if t, is known and the value of UAt is 
obtained by a trial and error substitution of 
several values of ¢, until the desired equality is 
obtained. 

By choosing six or more temperatures for t, 
and applying this procedure the point values of 
U At can be evaluated and graphical integration 
of equation (2) is possible, since values of UAt 
are known for increments of the total heat 


transferred. 
dq 
2 
| UAt (2) 


The method presented by Cotnuxn and 
HouGeN is accepted by all texts at the present 
time, but it is often pointed out that the method 
as applied is extremely tedious. There also 
appears to be no experimental evidence in the 
literature to support the design procedure. Smrru 
[i0}], however, does state that considerably higher 
heat transfer rates than expected were obtained 
on equipment designed by this method. 

Smiru attributes this to the method of cal- 
culating the heat transferred between points. 
This is made up of the heat of condensation and 
the heat transferred by the cooling of the gas, 
the uncondensed vapour and the condensate. To 
calculate the latter of these it is assumed that 
the condensate is at the same temperature as the 
main stream, which is not strictly true. Heat is 
removed by the cooling water as the condensate 
flows past the cooling surface and the condensate 
must be at a lower temperature than the mixture. 

Where the latent heat of vaporisation is small, 
SMITH suggests that the heat transferred by the 
cooling of the condensate is no longer negligible, 
as it is in the example worked out by CoLpuRN 
and Hovucen, and for organic vapour and gas 
mixtures some allowance must be made. 

SMITH gives a modification of the CoLBURN and 
Hovucen method whereby an estimate is made of 
the final condensate temperature from the 
empirical equation : 





ty — te! cop 
END 

From this ¢, is found and the total heat trans- 
ferred from inlet to outlet, is calculated, using ¢, 
as the condensate temperature, not t,. A plot of 
the terminal values of t, and ¢,, versus g is made 
and for a calculated value of t,, the estimated 
true ¢, is obtained by interpolation on this 
graph. The corrected heat transferred and 
corrected t, are then found and this new value 
of t,, used in equation (1) to solve for t, by trial 
and error as before. 

Sitver [9] gives approximate methods for 
different types of cooler condensers for estimating 
the gas film sensible heat transfer coefficient 
and the mean overall heat transfer coefficient. 

The calculations necessitate knowing the area 
of an existing condenser and its operating terminal 
conditions. Conversely the methods are applicable 
to the estimation of the surface area if the gas 
film sensible heat transfer coefficients are known. 

The methods are based on the assumption 
that h, is constant throughout the cooling process 
and that the ratio of the sensible heat transfer 
coefficient to the total gas film heat transfer 
coeflicient is equal to the ratio of the sensible 
heat capacity change to the total heat capacity 
change. 

However the calculation of the area required 
for a given set of conditions does not appear 
possible where the values of A, are not known 
for various equipment under different operating 
conditions, or where it is not possible to assume 
h, is constant. This may occur where there is 
a large velocity change of the gas-vapour mixture 
in the equipment caused by the removal of 
vapour from high concentrations of the vapour 
in the mixture. 

CoLBuRn [3] recently has presented an approxi- 
mate design method which involves only the 
terminal conditions of the condenser and elimi- 
nates the need for point to point heat balances 
and the trial and error calculation in order to 
obtain UAt. In this method the overall driving 
force from mixture to cooling water is expressed 
as an enthalpy difference (i, — i), and the 
conductance of heat is expressed as an overall 
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=) such 


coefficient per unit enthalpy difference | 
Cs 


that : 


dq/dA = (“/,.) i-< (4) 


CoLBuRN finds (‘ |) from equation (5) : 


WEL) = Wd) +1 Im) 


and he evaluates m, for the purposes of calculation 
at t.- 


' U : 
Assuming | ‘| ) constant, then equation (4) 
e 
becomes : 


(6) 


A = ——* _ 
(:/.) rm 


U , 
The value of ‘| however is known to 
c 
s 


vary and the arithmetic average of the top and 
U : 
bottom values of ( ‘| is used to evaluate A. 
c 
s 


Although the mechanism of the condensation of 
a vapour from a non-condensing gas is well under- 
stood, little appears to have been done to check 
experimentally the methods just described. It 
was felt that an experimental investigation on the 
validity of the design methods was most desirable. 
It is the purpose of this paper to present cxperi- 
mental data, obtained from a simple apparatus, 
to which design calculations can be applied and 
to compare the accuracy of the design methods 
with the known surface area used in the experi- 
mental work. 

Air mixtures with pure chloroform and with 
water vapour were chosen for the experimental 
work because of the differences in properties of 
these two vapours; notably the latent heat of 
vaporization. 


APPARATUS AND PROCEDURE 
The apparatus, as shown in Fig. 1 consists princi- 
pally of a boiler, a gas-vapour mixer, a verticle 
tube condenser, a condensate separator and an 
auxiliary condenser. 
The boiler is made from 1/16 in. steel sheet, 
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Fig. 1. 


10 in. 
The boiler is cadmium plated and is 


and the internal dimensions are 10 in. 

16 in. 
fitted with a gauge glass and safety lute. A } in. 
filling cock is provided and the exit vapour line 
leaves from an approximately 1 in. I.D. copper 
tube in the centre of the cover plate. A 3,000 watt 
Helicoil A.G.E., twin element immersion heater is 
sealed into the side of the boiler. The power 
input to the boiler is controlled by a 9 amps 280 
volt Variac on each element and is recorded on 
an Hvershed and Vignoles Recording Watt Meter. 

The vapour lines from the boiler are approxi- 
mately lin. I.D. and lead to the gas-vapour 
mixer which is also in glass. One foot before the 
mixer the air is introduced through a } in. glass 
line pointing in the direction of the vapour flow 
and the gas plus vapour enter the mixer tan- 
gentially. The mixer is constructed on the lines 
of a cyclone separator, the mixture of gas and 
vapour leaving from the top of the mixer. The 
mixed gas and vapour then pass to the condenser 
in approximately 1 in. glass lines; a three feet 
calming section is provided before the condenser. 

The condenser consists of a verticle double-pipe 
exchanger, with the gas-vapour mixture flowing 
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downwards inside the centre tube and cooling 
water flowing countercurrent to this in the 
annulus. The centre tube is of standard 12 gauge 
copper and is 1-289 in. O.D., 1-081 in. L.D. and 
the length of tube exposed to cooling water is six 
fect (or 1-70 sq. ft. internal area). 

The join between glass calming section and 
copper tube is a simple butt joint, both surfaces 
being first ground flat. The internal diameter of 
the glass at the join is the same as that of the 
copper tube. The two sections are held together 
by means of a bored rubber bung which fits into 
the annulus between the two metal tubes. The 
external tube forming the water jacket is standard 
8 gauge black iron pipe 2-055 in. I.D. 

A lin. LD., glass line which is six inches long 
leads from the bottom of the condenser and is 
fitted to the condenser in a similar manner to 
the calming section. The condensate formed is 
separated from the exit gas-vapour mixture in 
This is attached to 
the glass line which is flared at its end. Another 
glass line fitting under this takes off the exit 


the condensate separator. 


gas-vapour mixture. The condensate rate is 
dletermined directly from the volume collected 
in a known time. 

For the chloroform and air runs, an auxiliary 
condenser was used to remove as much chloroform 
vapour as possible before discharging to the 
atmosphere. A nine feet condenser with trichlor- 





ethylene cooled by dry-ice as the coolant, was used 
for this purpose. 

The maximum cooling water-rate available was 
pumped by a gear pump and metered through a 
rotameter. 

Standard thermometers are located at the top 
of the calming section, at the condensate outlet 
of the separator, in the exit gas-vapour line and 
in the inlet water line. Beckmann thermometers 
are used at the inlet and outlet of the cooling 
water to determine the water temperature rise. 
No correction was found necessary, at the water- 
rate used, to allow for fluid friction. 

For the system, chloroform and air, where the 
ratio of vapour to air is very high, the air was 
pumped by an Edwards Type IV Compressor (} hp 
1 kg/sq. em., 75 litres/min.) and measured 
through a calibrated constriction meter before 
entering the system. For the steam and air runs a 
Beecoxr Rotary Air Pump, 0-10 c.f.m., was used for 
the higher air rates needed in this case. The air was 
first filtered and metered through a previously 
checked } in. « } in. Venturi. For this system 
also the air was electrically heated before entering 
the vapour line. 

The boiler was calibrated for power input 
versus rate of evaporation for both water and 
chloroform. 
setting the equipment for total condensation 
without air flow and after allowing sullicient time 


This calibration was obtained by 


Table 1 





Vapour, Air Cooling 
Rate rate water 
System (Run| lbs lbs rate Inlet 
hour hour lbs 


1 47-2 5-80 2010 51:8 19-0 
Chloroform 2 47-2 5-80 2010 510 49-0 
and - 
air 3 47-2 5-80 2010 51-4 49-0 
4 47°2 5-80 2010 51:3 49-0 
| 1 1-1 | 2462 2010 748 | 763 
Steam 2 10-2 25-5 2010 750 75-7 
and - - : 
air 3 10-0 25:5 2010 74-6 a4 
4 wd 25-2 2010 750 75-7 


hour Exptl. | Caled. 


Gas- Vapour temperature “C Water Water | Con- | Con- 


inlet | temp. | densate | densate | Heat in 
Outlet temp. rise temp. lbs -- 
Cc Cc Cc hour Heat out 
Exptl. | Caled. 


36-3 34-4 13-7 1-14 30-0 32:3 0-04 


36-3 35-0 150 1-12 30-0 31-8 0-04 
36-4 b44 14-9 1-12 30-0 32:3 0-95 


| 36-3 34-4 14-7 1-11 30-0 32-4 0-96 


52°5 55-9 17-7 2-26 34-0 7-22) 0-08 
52-0 55-0 18-2 2-33 35-0 732) 0-06 
52 i) MO | 190 2-29 35-0 725) O97 
53-0 S44 191 226 | 345 727) «40-99 
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for the equilibrium the condensate rate was 
measured. A voltage regulator was found neces- 
sary to reduce the effects of mains voltage 
fluctuation. 

The boiler is lagged with a 1 in. asbestos sheeting 
and the condenser is lagged with 14 in. asbestos 
rope. All glass lines and the gas-vapour mixer is 
lagged with 3 in. of slag wool. 

The 


B.P. chloroform to remove the added alcohol and 


chloroform used was fractionated from 


then dried. The water used was taken directly 
from the mains. 

In making a run the water was turned on and 
the boiler set at the required vaporisation rate. 
Sullicient time was allowed for the vaporisation 
from the boiler to become steady and the air 
was introduced into the vapour line. The appara- 
tus was allowed to come to equilibrium and 
before readings were taken the mixer was drained. 
the 
a run was found to be negligible. 


mixer after 
The 


required to make a run was approximately two 


The amount of condensate in 


time 


hours. 

The results for both systems are listed in 
Table 1 together with the calculated saturation 
temperatures at inlet and outlet and the ratio 
of heat in to heat out, obtained by an enthalpy 
balance applied to the experimental results. 


APPLIED CALCULATIONS 


The design method of CoLBuRN and Hovcen [4] 
and also that of CoLBuRN [8] were applied to one 
set of experimental results from each system. Run 
numbers 4 in each case were chosen for the 
calculations, 

In applying the design methods the inlet 
temperature at saturation corresponding to the 
humidity of the inlet mixture was used for the 
purpose of the calculation and the outlet satura- 
tion temperature corresponding to the calculated 
humidity of the outlet mixture was also used. 
This was by the fact that the 
experimental temperatures did not correspond 
exactly to the temperature at saturation for the 
humidity involved, and that the Co_purn and 
HovuGen assumes cooling along the 
saturation line. In view of the fact that the 
actual temperatures did not differ greatly from 


necessitated 


method 





the theoretical temperatures as seen in Table 1, 
this procedure is justifiable. 

The calculations for both systems have been 
carried through simultaneously and for the sake 
of clarity sample calculations are given : 


Colburn and Hougen Method 
Table 2 lists the conditions known for each system 


Table 2 





Chloreform 


and air 


Steam 
and air 


0-869 
168-3 
129-9 


0-20 
120-2 
94-0 
Cooling water rate lbs /hour 2010 2010 
Inlet water temperature ‘( 14-7 19-1 
Pressure atm 1 1 
Inlet conditions : 

Vapour pressure atm 


Air rate lb moles /hour 
Inlet temperature °F 
Outlet temperature °F 


0-391 
0-609 


0-664 
Gas partial pressure atm 0-336 
Outlet conditions : 

Vapour pressure atm 

Gas partial pressure atm 
Ib moles vapour entering per hour 
Ib moles vapour leaving per hour 
Ib moles vapour condensed per 

hour 


0-151 
0-849 
0-558 
0-155 


0-382 
0-618 
0-395 
0-124 


0-403 





and the results of the overall mass balances. The 
inlet and outlet vapour pressures are obtained at 
the inlet and outlet temperatures respectively. 
Throughout the calculations the following 
properties were taken as constant : 
0-24 B.T.U, 
0-45 B.T.U, 
0-14 B.T.U. 
10 B.T.U. 
0-23 B.T.U. 


(ib) (-F) 
(Ib) (-F) 
(Ib) (CF) 
(Ib) (°F) 
(Ib) (“F) 


¢ for air 

¢ for water vapour 

c for chloroform vapour 
c for water 

c for chloroform liquid 
Molecular weight of air 29 
Molecular weight of chloroform 119-4 
Molecular weight of water 18 


The physical properties of the gas-vapour 
mixtures are given in Table 3 and were evaluated 
at the point temperatures chosen. The viscositics 
of the mixtures were calculated using the method 
of Bromiey and WILKE [1] while the molecular 
weights, densities and specific heats were cal- 


culated by standard methods. The diffusivities 
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| 

Viscosity 

Ib /(hr) 
(ft) 


Gas- Vapour 
temp. F 


Vapour | 
pressure 
alm 


System | Point 


: | 
120-2 


0-0305 


0-0310 


115-0 


Chioro- | 
form | 
and air 


110-0 


(4) 1035-0 0-0317 


(5) “100-0 0-0322 


(6) 


94-0 0-382 0-0329 


| i 
| Molecular | Density 


0-0295 


heat ft /hour 


0-356 





Specific | Diffusivity | Schmidt 
). 


Ne 


0-394 

| 0-439 
0-476 
0-521 
0-569 


0-631 


— -——-— 





168°3 0-0407 


160-0 


(1) 
(2) 


(3) 


(4) 
(5) | 


(6) 


0-391 


0-0419 


0-0430 


0-322 
O-254 


150-0 
0-197 


140-0 0-0439 
135-0 0-172 00-0443 


129-9 0-151 0-0445 


0-0539 0-322 


0-0563 


0-0588 


0-0612 


0-0624 


0-0634 


27°3 


| 1-14 
| 118 


1-12 


0-61 
0-62 


0-63 


0-63 


0-63 


0-63 
| 





were found from the equation of GrLLILanp [5] 
for chloroform and air but the value for the 
diffusivity reported by Suerwoop [8] for steam 


Table 4 


and air was used after correcting for temperature. 
The Schmidt numbers were calculated but since 
the variation for steam and air was slight an 





Tleat of 
condensation 
B.T.U./hour 


Points 


(1)-(2) 
(2)-(3) 
(3)-(4) 
(4)-(5) 
(56) 


Totals 


Chloroform and air 





2602 


2110 


Steam and air 


(445) | 


Totals | 7278 


Heat removed | Heat removed | 


from 
condensate 


from vapour 

B.T.U./hour 
26-1 
19:8 
15-9 
13-1 
12-4 


87-4 


27-8 

24-0 ; 

17:3 
73 


6-4 


Heat removed 
from air 
| B.T.U./hour 








Total 
B.T.U./hour 


1272 


2702 


2242 
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average value of 0-62 was chosen for calculating 
the mass transfer coeflicicnt. 

In Table 4 the results of the point to point 
heat balances are shown. 


e.g. Points (4)-(5) Steam and Air 


Heat of condensation 
0-197 0-172 


0-869 _ - 
* 10.803 ~ 0-828 


) 18x 1014 = 584 B.T.U., 


Ifeat removed from condensate (con- 
densed assumed to leave at same 
temperature as the mixture) 

0-391 0-172 

0-609 0-828 


0-869 x ) 18 x5 33-9 B.T.U., 


Heat removed fron uncondensed vapour 
0-172 


0-828 x 86918 xk O45 x 5 


7-3 B.T.U. hr 


Heat removed fron air 
= 0869 x 20 x 024 x 5 


= 30-2 B.T.U./hr 


Heat transferred between 


(4) and (5) 


points 


655 B.T.U./hr 


The total heat transferred between inlet and 
outlet is thus 3781 B.T.U./hr for the chloroform 
and air system and 7795 B.T.U. /hr for steam and 
air. This enables the water temperature at each 
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point and the total water temperature rise to be 
calculated. Table 5 gives these values. 
The mass velocities and Reynolds numbers are 
Table 5. 





| Total heat! W 
transferred bo 
B.T.U./ ¢ 


Water 
temp. 
© F 


System Points 





At (1) 


(1)-(2) 


(2)-(3) 


(3)-(4) 
(4)-(5) | 497 


(5)-(6) | 499 


Totals 3782 





Chloro- 
form 
and air 








0 


| 


1-34 
1-12 
2 


0-82 


~(1)-(2) 
(2)-(3) 


(3)-(4) | 1646 — | 


2702 

2242 
Steam 
and air 


(4) (3) ae 0-33 
(4)-(6) 
Totals 


0-27 


3°88 | 








System 


sane ancora 


25,400 


8,310 


5,370 
Chloroform and air 


19,300 | 


j 
0-00382 


0-0033 






































Steam and air 








4,470 





(6) 4,390 
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shown in Table 6 for each point. These are 


evaluated as follows : 
e.g. Point (2). Chloroform and Air 


1-081 in. 
0-00637 sq. ft. 


Tube diameter 
Cross-sectional arca 


0 
p= O20(119-4 » 


0-00637 
6,530 Ib /(hr) (sq ft) 
0-0305 Ib /(hr) (ft) 
= 1081 «x 6530 
“y2 x 0-0305 


aid pg at (2) 
So N 
Re 19300 

The j factor is read directly from figure (1) of 
reference [2]. 

The values of the sensible heat transfer coelli- 
cient, h,, and the mass transfer coefficient, K, are 
calculated from the equations of Cuiuron and 
CoLaurn [2]. 


(7) 
jG 


4 
aM. Pes - p} 


value of 0-80 is assumed for the 


K = 
The Prandtl 
number for chloroform and air giving 

(Pr) 


For steam and air a value of 0-76 is taken for the 


0-862 


Prandtl! number giving 
(Pr) 
The calculated values of A, and (A 
are also listed in Table 6. 


0.833 


Pgs) 


e.g. Point (2). 
j = 0-0033 


Chloroform and air. 
G = 65380 
cp\ 4 . 
c = 0-180 (F) = 0-862 
From equation (7), 
00-0033 =« 0-180 x 6530 
~ —_ 
= +50 B.T.U. /(hr) (sq ft) CF) 
Also M,, = 83-2 
“ 


h 


= 0-439 
pb 
From equation (8), 
0-0033 x 6530 
B32 X Per x (0-439)8 
0-448 
Pat 


K = 
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Heat Transfer Coefficients other than the Gas Film 


The heat that is transferred to the surface of the 
condensate film must pass through the resistances 
of the condensate film, the tube wall, scale and the 
water film before being removed by the cooling 
‘yr nd 

These resistances calculated as 


water. are 


follows : 


(a) Condensate film, h, 
Values of 4, are calculated from the equation of 
Nussecr [6], 


(9) 


h, = 0.98 k ote 
aT 


For chloroform and air h, was calculated at each 
point. For steam and air, where the condensing 
film coefficient is much higher, the value at the 
bottom of the condenser, where the resistance is 
greatest, was taken as constant throughout the 
condenser. This value was calculated from 
equation (9) as 1360 B.T.U./(hr) (sq ft) (°F). 
Table 6 also lists the point values of h, for 
chloroform. 


(b) Tube Wall, h, 
220 « 12 


m ~~ 0-104 
> 25,400 B.T.U./(hr) (sq ft) (°F) 


h 


(c) Seale, hy 


A scale coeflicient of 1,000 is allowed for 


both systems. 
hy = 1,000 B.T.U./(hr) (sq ft) (°F) 
(d) Water Film, h, 
The film coefficient the water side was 
evaluated from CoLBuRN’s equation for fluids 


inside tubes in turbulent flow, recommended by 
Perry [7] for turbulent flow in an annulus : 


(%) (poem (40) 


The value obtained is 


h,, = 179 B.T.U./(hr) (sq ft) (°F) 


on 


(10) 


When these coeflicients are corrected to the 
inside area of the condenser tube, 


h, = 179 B.T.U./(hr) (sq ft) (°F) 
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for chloform and air and In the same way values of U At, U and h,, have 
h, = 159 B.T.U./(hr) (sq ft) (°F) been calculated for the other points. These final 


o 
-~ 


values are listed in Table 7 together with the 


for the steam and air calculation, since a constant 
values of q and h,. 


value of h, 1,360 was used for steam. 





10 
Evaluation of U At 
The values of UAt at each point were calculated ; 
by substituting in equation (1) the previously 


determined variables, and solving for t, by trial — - — 
and error. 


- 4 





e.g. Point (2). Chloroform and Air 








4 
| 

— 

311 


h, 114 B.T.U./(hr) (sq ft) CF) ‘ =e 

For ¢, 87-3 F equation (1) becomes, 7 TL OROFORM- Aah 
0-445 

4:5(1150 87-3) 4 119-4 & 108-5 (0-600 —O0-330) 
( 3 


-° 





52! 
114(87°3 50-7) U (15-0 59-7) 
$121 3146 = 553 U 
Average U At = 3130 B.T.U./(hr) (sq ft) 
U 56-6 B.T.U./(hr) (sq ft) CF) 
and the hypothetical coefficient is 








4 i i 


3130 1000 2000 3000 4000 $000 
- o.< q 8ru/m 
(115 — 87-3) 


113 B.T.U. (hr) (sq ft) (OF) Fig. 2. 











System Point 


5030 
O-86 3130 56-6 d 2509 
Chloroform b 57 0-683 2400 17 1658 
and air 
0-621 0-565 1830 7 9906 


0-660 O-482 1420 34- 52-3 100 


0-704 0-393 { 1000 


0-737 | g120 
791 37 M 6260 
Steam and : 0-842 


air 
O-877 


O-sol 1 


0-906 0-08 
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From a plot of 1/U At versus q, as in Fig. 2, 
equation (2) is graphically integrated to give the 
surface area required. The area obtained for the 
chloroform air system is 1-68 sq ft, which is 
equivalent to a condenser tube length of 5-94 ft. 
For the steam and air system, the area is 1-63 sq. 
ft., or 5-76 ft tube length. 

The good agreement between the calculated 
areas and the expected area of 1-70 sq ft is 
apparent. It suggests that any further modifica- 
tion, such as that of Smirn is unnecessary, for 
the calculated area for the chloroform and air 
run agrees to within 2°%, of the expected and the 
calculated area for the steam and air run agrees 
to within 5%. Furthermore the modification of 
SMITH would lower the calculated area which 
would increase the deviation between calculated 
and expected areas. For these reasons the Smiru 
modification was not applied, nor was any other 
allowance made in the calculations for the cooling 
of the condensate. 

The average condensate temperature, ¢,, at the 
bottom of the condenser, obtained by taking the 
arithmetic mean of the wall temperature, ¢,,, and 
the condensate surface temperature, f,, is calculate 
as 71-5°F for the chloroform and air system and 
80-3°F for the steam and air system. These 
temperatures are lower than the expected 
average condensate temperatures of 86-0°F 
(30-0°C) and 94-1°F (34-5°C) for the respective 
chloroform and water condensate. This is 
perhaps due to reheating of the condensate 
by the gas-vapour mixture before it can be 
separated from the mixture in the experimental 
apparatus. 


Colburn Approximate Method 


The approximate method of CoLtpurn [3] out- 
lined previously was also applied to both systems. 
The results of the calculations for this method 
are given in Table 8. For the value of h, used 
here the effect of "/h, was quite considerable 
and in this case the accuracy of the method 
depends in part on the exact evaluation of the 
slope of the enthalpy versus temperature line. 
Actually the most difficult part of this calculation 
is the determination of m. The value of m in 
Table 8 was found as follows : 





1 
A plot of log p, versus 7 was made, as 


suggested by the Clausius-Clapeyron equation, 
over the temperature range involved in each case, 
i.e. from inlet to outlet ¢,. From the resulting 
dp 
Po The value of 
dt 


straight line was found. 


Po 
is also known from the humidity-vapour pressure 
equation and hence by application of the function 


; du - 
of a function rule results. Since m is evaluated 


at t,, as suggested by CoLburN, m becomes 
Differentiation of equation (11) 
i = Agy H + ¢, (t — 82) (11) 
with respect to ¢. remembering that 


c=c Le <H (12) 


s (attr) (vapour) 


oo dH . , 
and substitution of : gives an equation for the 


l we 
value of “" between the temperature limits 


chosen. 
This method for evaluating the slope of the 
enthalpy line was found to be much simpler to 





Table 8 
| Chloroform Steam 
and air 


and air 
>< | — 


Top | Bottom Top | Bottom 


Vapour temp., ty . | 120-2 | 940 |168-3 | 129-9 
Water temp.., ¢,, | 60-3 58:5 | 703 | 66-4 
Vapour pressure, — 0-664 0-391) 0-151 
Enthalpy, i, B.T.U./lb 1070 | 333 485 147 








232 «#(458 «| 128 





i,,* jee 101 26-6 | 23-5 
Ai |96 | 
m at ty | 205 570 | 793) 1-06 
h, | 587) 802) 813) 5-74 
t (188 B 0-60 F bt | O29 
h, (179 j1s0 159 
(Us/e,) | 2-37 | $89 | 985) 150 
———__ — — —_—_— | ae 
Total B.T.U./hr 3781 | 7795 
(U,/c,) arith. average Bs | 12-9 
(Ai)s mm. 515 255 
Calculated area sq ft 2-35 2-40 
Equivalent tube length ft 8-3 85 

| 








186 
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use than standard mathematical methods, such 
as the Differentiated Gregory- Newton method and 
the Legrange Interpolation Method. As long as 
the temperature range involved is relatively small, 
so that the variation in latent heat of vapourisa- 
tion is small, then the method is quite applicable. 
Using this method the following equation was 
obtained : 
di =K'p,M, f{ 


dt (1—p,)* 7? M, | 


where for chloroform and air K’ = 8030 


Ase + ¢(rapour) * (¢ — 32)} + Ce 


and for steam and air K’ = 8990 


The final value of the area required as shown 
in Table 8 gives areas for both calculations 
approximately 40°%, greater than is needed, which 
is reasonable considering firstly the assumptions 
on which the method is based and secondly the 
fact that an oversize area is obtained. For the 
cases studied the method is certainly conservative 
from the point of view of capacity. 


CONCLUSION 
The design procedure put forward by CoLBURN 
and HouGeEN applied to known conditions has 
been found to give an area which is in good 
agreement with the expected area for the systems 
chloroform and air, and steam and air. The 
modification of this method proposed by Smiru 

has been found unnecessary. 
The approximate method of CoLBuRN applied 
to the same conditions has been shown to give an 
area which is conservative from the point of view 


of capacity. 
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NOTATION 
A = area, sq ft 
¢ = specific heat at constant pressure B.T.U./(Ib) (°F) 
¢, = humid heat B.T.U./(Ib) (°F) 
d = diameter of tube, ft 
d, = equivalent diameter, ft 
D = diffusion coefficient, ft®/hr 





dA = element of surface area 
dq = increment of total heat transferred per unit time 
£1, = acceleration due to gravity 4:18 x 10° ft/(hr)* 
G = mass velocity (Ib/(hr) (sq ft) 
h = film coefficient of heat transfer, B.T.U./(hr) (sq ft) 
(*F) 
hg = combined water, metal wall and scale con- 
ductances 
= condensate 
dirt or scale 
metal wall 
combined conductances other than the gas 
film 
= gas film 
hy = dg/(tg — te) dA 
= water film 
= humidity of mixture, Ib/lb dry gas 
i=enthalpy of saturated gas-vapour 
B.T.U./lb of dry gas 
= enthalpy of saturated gas-vapour mixture at (, 
B.T.U./lb of dry gas 
= enthalpy of saturated gas-vapour mixture at the 
temperature of the cooling water B.T.U./Ib dry 
gas 
= logarithmic mean of i, and i,,* 
j = heat transfer or mass transfer factor 
= thermal conductivity, B.T.U./(hr) (sq ft) (°F /ft) 
. = molar mass transfer coefficient, Ib moles/(hr) 
(sq ft) (atm) 
= length, ft 
= slope of curve of i,,* versus f 
= molecular weight 
M, = air 
M,, = average gas-vapour 
M,, = condensable vapour 
M, = inert gas 
= partial pressure, atm 
Pe = Vapour pressure at (, 
Pg = non-condensable gas partial pressure in the 
main body 
Pg’ = that adjacent to condensate surface 
Pgs = logarithmic mean of p, and p,’ 
P, = vapour pressure 
q = heat transferred, B.T.U./hr 
T = absolute temperature, °R 
t = temperature, °F 
= average condensate temperature 
= condensate surface temperature 
== gas-vapour mixture temperature 
= wall temperature 
= water temperature 
At = overall temperature drop 
U = over-all coeflicient of neat transfer, B.T.U./(hr) 
(sq ft) (°F) 
U,/¢,= total gas-film heat transfer coefficient per unit 
enthalpy difference, B.T.U./(hr) (sq ft) (B.T.U./ 
Ib) 


mixture, 
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I’ = mass rate of flow, Ib/(hr) (ft of wetted perimeter) REFERENCES 


{1] Bromuiey, L. A. and Wikre, C. R.; Ind, Eng. Chem. 


A = latent heat of vapourisation, B.T.U./Ib 
1951 43 1641-8. [2] Cuitton, T. H. and Coisurn, A. P.; 


pp = viscosity, Ib, (hr) (ft) Ind. Eng. Chem, 1934 26 1183-7. [3] Cotpurn, A. P. and 

Inst. Mech. Engrs. (London) Proc. 1951 164 448-58. [4] 

p = density, Ib/cu, ft. Coisurn, A, P. and Houcen, O. A.; Ind. Eng. Chem. 

ee 1934 26 1178-82. [5] Perry, J. H.; Chemical Engineers’ 

Pr k’ Prandtl number Handbook, 3rd Ed., p. 538, New York, McGraw-Hill Book 


Co, 1950, [6] Perry, J. H.; Ibid. p. 476. [7] Perry, J. 
. » . H.; Ibid. p. 470. [8] Suerwoon, T. K.; Absorption and 
So. ae oe Extraction, Ist Ed., p. 21, New York, McGraw Hill Book 
Co. 1937. [9] L. Sitver; Inst. Chem. Engrs. (London) 
Trans. 1947 25 30-42. [10] Smrrn, J. C. ; Ind. Eng. Chem. 
34 1248-52 (1942). 


Nr , Reynolds number 
rv 
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Summary —A graphical statistical method is developed for the study of sulphur compounds in 


mineral oils. 


The corresponding graphs are based on the ultimate analysis and physical constants and 


2 
: : , , - , n* — 1 
functions such as refractive index (n), density (d) and specific refraction ( z ° ). 
n 2 


A statistical graphical representation of the free enthalpy of formation of n-thiols, n-alkenes 


and H,S, n-alkanes and sulphur is given. 


Résumé 


binaisons sulfurées, présentes dans les huiles minérales. 


L’auteur propose une methode graphique et statistique pour la recherche des com- 


Les figures correspondantes sont basées 


sur l'analyse élémentaire et sur des constantes ou fonctions physiques comme n, d et 


n?—1 1 
n242 4° 


L’auteur propose également une réprésentation graphique et statistique pour le potentiel 
thermodynamique moléculaire de formation des n-thiols, alkenes et H,S, n-alkanes et soufre. 


INTRODUCTION 


To arrive at a better understanding of the 
composition of mineral oils the study of the 
sulphur compounds occurring in practically all 
these oils is of importance. 

A preliminary investigation is made into the 
physical and chemical properties of sulphur 
compounds and of their mixtures either in the 
presence or in the absence of hydrocarbons. 

A new method is proposed for improving our 
knowledge of sulphur compounds occurring 
in mineral oils and a thermodynamic study is 
made concerning possible conversions of sulphur 
compounds, 


SECTION 1 
About a year ago D. W. van Krevecen [1] 
applied a new method in coal research by using 
graphs, based on the ultimate analysis of coal. 
Such graphs can be also applied to mineral 
oils. 

The majority of the mineral oils contain 
considerable quantities of sulphur compounds 
besides hydrocarbons. Compounds containing 
other elements such as oxygen and nitrogen in 
addition to carbon and hydrogen are generally 





only present in smaller quantities in mineral 
oils. 

It is important to know which group of sulphur 
compounds is dominating or preponderant in the 
oil as one type may behave quite differently 
from the other. The type of sulphur compounds 
will have also a specific influence on the properties 
of the crude oil. 

In Figs. 1 and 2, graphs are given for compounds 
which consist of C, H and S. These graphs can 
be applied to pure sulphur compounds or to 
mixtures of sulphur compounds, They enable 
the contruction of lines which represent the 
different groups of sulphur compounds. 

Each defined sulphur compound is represented 
in the graph by a point, e.g. thiophene (C,H,S) 
in Fig. 1 by H:C=1-00 and 8:C = 0-25 
and in Fig. 2 by C: H = 1-00 and S:H = 0-25. 

In Fig. 1 hydrogen sulphide (H:C = © 
and §:C = ©) lies in the infinite. 

Fig. 2 has the advantage that H,S(C : H = 0.00 
and S: H = 0-50) is represented by a point and 
that hydrogenation and dehydrogenation can 
be represented more clearly. 

Two examples might illustrate how to use the 
graphs (c.f. Fig. 2) : 
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From H,S (C:H 40 





= 0-00, S:H= 0-50) 





and C,H, (C:H = 





0-40 and S : H = 0-00) 





thiophene (C: H = 





3.5 
1-00 and S:H = 0-25) 





can be obtained re- 








sulting in a dehydro- 
genation. The equi- 





molecular mixture of 9396 








hydrogen sulphide 





and butane is repre- 











sented by the point 





C:H = 0-33 and 





S:H = 0-088. From 25 











H,S (C : H = 0.00 and 











S:H = 0-50) and 











alkenes C,H,, (C:H 45 





= 0-50 and S:H = 





0-00) thiols and 














thioethers can be 








hydrogen sulphide Z 
| 


formed, which are re- 























presented by points 15 |_| 
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such a mixture. By 

































































determining the re- 
fractive index, density 
and molecular weight 
(n — d — M method) even an accurate analysis 
should be possible. In such an investigation 
of crude oils the graphs discussed above could 
be of good help. 

Suppose that a crude oil fraction containing 
sulphur compounds is exposed to an operation 
so that the composition is changed as a result 





1.0 15 2.0 25 
Fig. 1. 


of this. The investigation of the original fraction 
could be done as follows : 

The fraction should be separated in two parts, 
i.e. hydrocarbons and a mixture of sulphur 
compounds. The hydrocarbons can then be 
investigated by means of the n-d-M method, the 
sulphur compounds by means of the graphs plus 
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the determination of the refractive index, density, 
molecular weight and possibly other physical 
constants, the sulphur content and the ultimate 
analysis. 

The operation with the original oil fraction 
is then carried out and the resulting product is 
investigated in an analogous way. 

Thus it will be possible to obtain a certain 
knowledge of the behaviour of the hydrocarbons 
and of the sulphur compounds during the process. 
It might be possible that during the process 
sulphur compounds have been transformed into 
hydrocarbons and other sulphur compounds, 
whilst also hydrocarbons might react with 
sulphur compounds thus giving other hydrocar- 
bons and sulphur compounds. 

The main difficulty is the complete separation 
of hydrocarbons from sulphur compounds. Those 
sulphur compounds containing a high percentage 
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of carbon and of hydrogen resemble very closely 
physically and chemically the corresponding 
hydrocarbons. There are indications that, in 
certain cases, a complete separation of hydrocar- 
bons and sulphur compounds seems _ possible. 

As a secondary problem there remains in certain 
instances the presence of oxygen and nitrogen 
compounds. 

In the processes of synthesis of sulphur com- 
pounds from hydrocarbons and hydrogen sul- 
phide, or from hydrocarbons and sulphur, or 
from hydrocarbons and sulphur compounds, 
the mixtures formed will not have such a 
complicated composition. Therefore such a 
synthesis can be more easily studied. 

As the physical constants of a large number of 
sulphur compounds are already available it is 
possible to prepare a graphical representation 
as shown in Fig. 3. In this graph the specific 
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Fig. 3. 


refraction ” - is plotted versus the percent- 
n?+2 d 
age of weight of sulphur. 

On the lines the values for n7Z are indicated 
which enables an easy orientation. Fig. 3 is 
only of value if special types of compounds, e.g. 
strongly unsaturated aliphatic thiols have not 
to be taken into consideration. 

The accuracy in the application cannot be 
compared with the n-d-M method for the analysis 
of hydrocarbon mixtures. 

Nowadays an intensive study of 
the 
takes place all over the world. 
question can be therefore improved gradually. 

The present knowledge of the physical con- 
stants of many 
enables us to predict by inter-or extrapolation 
physical constants (x7), d%, boiling point) 


sulphur compounds not yet investigated [2]. 


sulphur 


compounds sponsored by oil companies 


The graphs in 


sulphur compounds already 





SEcTION 2 


For a more detailed study of the behaviour of 
sulphur compounds in mineral oils in the chemical 
and physical processes thermodynamical con- 
siderations cannot be omitted. 

With the aid of the knowledge from the 
literature 


molecular 


[3] it is possible to calculate the 
free enthalpy of formation per C 
atom (AG, in kilocalories per mol. per C atom) 
from thiols, (graphs in Fig. 4). 

The molecular free enthalpy of formation is 
plotted versus the number of C atoms. The 
n-thiols have been compared with the n-alkenes-1 
and hydrogen sulphide and with the n-alkanes 
and sulphur. The graphs are valid for the 
gaseous phase at one atmosphere. 

For the temperatures 300°K, 500°K and 
800°K, respectively, the curves are calculated 
and given in Fig. 4. From the place of the curves 
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in the diagram the stability of the considered 
compounds can be derived. 

At 300°K the normal thiols are stable compared 
with the corresponding n-alkenes-1 + hydrogen 
sulphide and stable compared with the n-alkanes 
+ sulphur. 

At 500°K the thiols are just metastable com- 
pared with the systems mentioned. If the rate 
of reaction would be high enough the normal 
thiols could decompose in n-alkenes-1 + hydrogen 
sulphide at 500°K and 1 atmosphere. Thus 
could be explained the formation of hydrogen 
sulphide in the final distillation of refined crude 
oil fractions which still contain sulphur 
compounds. 


J. An oe ae 








From the graphs it can also be concluded that 
alkanes can be transformed by dehydrogenation 
with sulphur into alkenes. 

These graphs have not a quantitative value 
as their basis is not quite certain. In the near 
future it may, however, be expected that more 
thermodynamical data will be available, as 
recent literature shows thermodynamical pro- 
perties even of higher sulphur compounds in the 
preparation of very pure sulphur compounds. 
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Book review 


Book review 


Dr. Inc. Wine tm Guz: Vergasung fester Brennstoffe, 
Springer-Verlag, Berlin, Gottingen, Heidelberg, 1952. 102 
pp. 4 illustrations. IDM 9.60, 

The author has typified the matter dealt with in his book 
by the sub-title, Stoffbilanz und Gleichgewicht - Eine 
Darstellung praktischer Berechnungsverfahren. 

The author restricts himself mainly to a discussion of 
the various methods developed for calculating the equili- 
brium composition of gases that may be produced when 
solid fuels are converted with gasifying agents. 

Moreover, these calculation methods are illustrated by 
examples in which even the formation of sulphur-contain- 
ing compounds, such as H,S, COS, CS, has been reckoned 
with. It is regrettable that no use is made in these examples 
of the heat balances that may be drawn up, although the 
drawing up of such balances in general is dealt with in a 
separate section. 

It is owing to the omission of the heat balances, that 
the results of some of the calculation examples given differ 
considerably from those obtained in practice. 

This applies for example to the first (and very simple) 
example of a process according to which pure carbon is 
gasified in counter-current with dry air at atmospheric 
pressure and the equilibrium composition of the gas is 
calculated at a temperature of 650°C (p. 70). However, in 
adiabatic systems the equilibrium temperature in the 
gasification of coke, carried out with dry air of room 
temperature, is easily twice as high. 

This explains why in practice it is quite possible by 
means of cold dry air to produce a gas containing 
less than 1 per cent of CO,, instead of the 9-89 per cent 
stated by the author. 

Those who should feel inclined to adopt the calculation 
examples given here as a guide for carrying out equilibrium 
calculations, must be well aware that in order to obtain 
practically useful results, the character of the gasifying 
agent on the one hand and the height of the equilibrium 





temperature on the other hand may not be selected 
independently of each other, but that the relationship 
between these two factors must be calculated with the 
help of a heat balance. 

As the author has prefaced his main subject by a rather 
extensive general introduction to the gasification theory, 
his bcok deserves more general interest than might be 
aroused by the sub-title. The introduction is certainly 
worth reading, although the critical reader will not be 
entirely satisfied by it. 

For example, many will rightly have doubts as to the 
truth of the statement that in practice the Boudouard 
equilibrium and the water gas equilibrium should always 
be attained (p. 22), although it must be admitted that the 
value of the equilibrium calculation for the gas technician 
can be hardly overestimated. 

In the chapter on the kinetics of the gasification reaction, 
the author emphasizes the great difference between the 
carbon temperature and the gas temperature that may 
occur during the gasification process, as well as the influence 
of the physical transport resistance on the rate of the 
gasification reaction. 

However, the remark made in this connection to the 
effect that at temperatures of from 700-750°C the physical 
resistance should already far exceed the chemical resist- 
ance is certainly not acceptable. Generally speaking, this 
only holds for the reaction between carbon and oxygen 
but not for reactions between carbon and carbon dioxide 
or steam. 

Although some more remarks of this kind might be 
made, this does not alter the fact that it is very desirable 
to possess a book in which the various methods for calcu- 
lating the gas equilibria are collected and which contains 
an introduction in which an expert confronts the reader 
with the more fundamental problems connected with the 
gasifying process. 

W. van Loon, 
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